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Cl .T y no. or g 
* ^ Elements 0 


no. 

of 


01 juegrees 
freedom 


Roll force 
(tons) 


1 

80 
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2 

2k0 

310 

3 

320 

iflO 

if 

if 00 

510 


39.^215 

23.9^6 

22.6285 

18.9730 


Experimental Value of RoH Force = 16*56 tons. 

It can be seen that the present results with 
^00 elements are very close to experimental values . The 
discripency between experimental and theoretical values 
can he reduced by taking more number of elements and by 
taking elasticity of the rolls into account ^ Fig. I^f). 
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ABSTRACT 


In this work the finite element laethod has 
"been applied for the analysis of cold rolling process. 

The prohlem is solved by dividing the material into 
plane strain triangular elements. Two models of elastic- 
plastic hehavior of the material, namely, the elastic- 
perfectly plastic and the Ramterg Osgood models, were 
considered in this work. The analysis provides the roll 
forces, roll torque, pressure distiihution on the work- 
roll interface a.nd the deformation pattern and flow of 
material under the rolls. The results obtained were 
found to be in good agieement with the available experi- 
mental results . 

Further, the mathematical programming tech- 
niques >;ere applied to design the cold rolling process 
optimally. The design problem is formulated using 
homogeneous theory of metal deformation due to Eland 
and Ford, The follomng three optimisation problems 
were solved; 

(i) the design of single pass rolling process 

■i'jith the deflection of rolls as the objective 
for minimisation, 

(ii) the design of single pass rolling process 

with power as the objective for minimisation, 



(iii) the design of nuLtipass rolling process 

T'jith the sum of deflection of rolls in all 
the nills as the objective for minimi sat ion. 

It is found that tensions play a more impor- 
tant role in tiio case of rolling of thinner strips 
compared to that of tliicker strips . Further, suppl^/ing 
power through the rolls has been found to be a more 
efficient way of transfering the poirer for the defor- 
mation of metals. In the case of a tandem mill, most 
of the reduction has been observed to talie place in the 
first mill itself. 



:i-I-ZEODUCTIOIT 


1.1 NATURE OF FETAL DEFORI'iATI OK PROCESSIKG PROBLEM 

The mechanLca.1 processes of greatest interest 
from the standpoints of total tonnage handled and 
technological advances made, are those involving steady- 
flow, that is, where stresses, strains and velocities 
at some particifLar reference location do not change 
with time. Examples include rolling^, drawing and 
extrusion* The finite elecBnt analysis of rolling pro- 
cess and the design of cold rolling mills is considered 
in this work. 

As Coffin [^lj in his re-view puts it, ”the 
"Ultimate obiecti-ve in mechanical deformation process- 
ing is. the attainment of a unified theory which would 
predict all important paranBters designated by -the 
equipn^nt supplier, the producer and the consumer of 
the product," On the basis of initial properties and 
controllable variable^- the theory should pro-vide the 
force, pressure distribution and power requiren^nts 
necessary for the construction of equipment, Next it 
should describe point to point stress, strain and 
strain rate relationships during the deformation to 



enable the producer to predict the residual stresses, 
the final shape of the product and whether failure 
would occur. Thirdly and finally, it should provide 
the necessary’- information leading to a prediction of 
the final properties of the product, including mecha- 
nical properties, texture and anisotropy. 

The proposed r;;athem3.tical model of rolling 
process, namely, the finite element model provides 
the forces, pressure distribution and torque on the 
basis of initial properties and va,richles and des- 
cribes the stresses and strains at various points . 
Finally'’ bj^ knoimLng the strain histor 3 f and stresses, 
metallurgists will be in a position to pr-edict the 
final properties of the product. 

Once the forces, pressure distribution and 
torques are known, it is the designer's job to design 
the machine with the main objective of Eiinlmising the 
cost (initia-l as well as running costs) and / or of 
improving the quality of the product. In metal defor- 
mation processes where large production rates are 
involved, the running cost assumes more importance 
than the initial cost. In most cases, where auto- 
mation is there, the cost of running is almost 
directly proportional to the power i«quirements of 
the machine. For the improvement of quality, it is 



desired that the machine should v^ork near the design 
configuration, i.e^, the deformation of components 
should he small. In general, the minimisation of the 
poorer and the minimisation of the deformation can not 
he achie-ved simultaneously. In fact one can he 
achieved only at the cost of the other. 

In practice, there vdll he certain restric;«!i 
tions on the performance of the machine. These res- 
trictions limit the choice of the design parameters. 
Hence the designer has to choose the design para- 
meters so as to .sati''fy il-l th<= frictions placed 
on the performance of the machine, which also minimise 
the po^i?er requirements of the machine or the defory 
mation of the components. In the present work, the 
mathematical prograrmiing techniques are applied to 
find the design parameters which will achieve 
objective stated earlier under a specified constraint 
set for cold rolling of strips , 

1»2 LI2ERATUHS SURVEY 

1*2.1 Review of Cold Rolling Theories 

- ;ioiiujnr4jixx."- 

Tl^ homogeneous flow theory in cold rolling 
assumes that th$ stress, strain and velocity at all 
points on a plane normal to the axis of the bar are 
constant. This assumption, along with a suitable 



yield criterion, can be used to derive the differen-^ 
tial equation for the equilibrium of a sgall transverse 
slice of the total cross section of the metal being 
deformed, shovn in Fig, (1) (valid for wide bars, • 
'sheets and strips), as . 

d (fh) = sdh+2Asdx = 0 (1.1) 

This equation was first developed by Von- 
Eaji3an{^2j , Depending on the assumed degree of 
simplification of the geometry, nature of filctional 
effecis and material properties, it is possible to 
solve the differential equation and obtain solution 
for the pressure distribution along the roll -work 
interface, the total roll force and the torque re- 
quired. Important theories that have evolved froBi 
this approach are those of Von-Karman £ 2 ’] , Trinks 
£3], Tselikov 3 5 Nadai £5] . and . Bland and Ford 

In this study, Bland and Ford£6j theory is used • 
for the optimization of the rolling process, as it 
is found to give good results -vdien compared with the 
e.xperimental results £21 J, These theories predict 
■ the important paiameters needed for the construction 
of equipment but, because of the assumption of homo- • 
geiasous flow, are of no value in predicting the stress 
strain distributions during or after the defer- 



The experimental investigations made by 
I'iac Greger and Coffin £7} » Hundy and Singer 
show that non-homogeneity is present in rolling 
process and it is important to take it into account. 

The slip line -field approach is based on the very 
feature which produces non-homo geneous deformation 
cind, hence ^ is a powerful analytical tool in metal 
forming processes, Eut according to Hill 9 j accu- 
rate slip-line solution of the problem of strip 
rolling is not yet known and even the qualitative 
appearance of the slip-line field has still to be 
demonstrated. It is assumed to be applicable to 
situations where reductions are small to permit one 
to assume that the arc of contact is a straight line. 
Further, , the appearance of the slip-line field is 
ass-umed to correspond to the drawing of sheets 
through a smooth wedge-shaped die, which gives a 
moderate reduction in thickness 

Both the homogeneous flow and the slip-line 
theories are applicable to plane-strain problems only. 
When it comes to the analysis of rolling of rectangular' 
bars, where the spread cannot be ne^ected, the theo- 
3?ies have hardly progressed beyond the empirical sta'gfe* 




1*2.2 Be view of KLnite Element I^^thods in 
Metal Forming Processes 

The finite eleinent method vras developed -for 
the analysis of complex structures and is in the 
process of rapid extension to a. wide range of non- 
structural' prohlems . Metal forming can he considered 
as a non -structural problem*. The finite element' 
method has been applied to simple cases of metal 
forming like upsetting [loj and indentation QnJ » 
where studies were made on the development of the 
plastic zone, the load-displacement relationships 
and the stress and strain distributions . More re- 
cently finite element method is applied to the , 
problems of plane -strain and axLsyrsnetric extrusion 
£ 12 ^ where the effect of frictional coefficient on 
the spread of plastic zone, the pressure displace- 
ment curve and the stress and strain distributions 
are studied for the initial non-steady state extru- 
sion using incremental -load method of aralysis . 

In this study, the cold rolling of strips 
is analysed as a plane strain problem using the 
fihite element method. The me'thod is equally appli- 
cable to problems where spread can not be neglected. 
The incremental displacement iiethod is used for the 
non-linear analysis of this problem. 




1 .2.3 Revlev of Design of Rolling Processes 


The ^problem of design of rolling mills has 
always been of interest to mechanical engineers for 
a very long time. This is evident from the fact 
that many theories and empirical relations are avai- 
lable to calculate roll forces, pressure distribution 
and power requirements of the deformation process in 
cold rolling |^2 - 6 J . These quantities are very 
important for the design of rolling mills and / or 
processes. But nothing much has been done in the 
direction of finding optimal conditions and parame- 
ters for the rolling process. Avitzur’s article 
(jkj is one of the few attempts made in applying some 
theoretical considerations to the problem of tandem 
mill optimization in hot rolling. In this article, 
the rate of production is maximised. In this problem, 
the inter stand tensions, the thicknesses of the 
material in between any two adjacent mills and the 
speeds of the mills are taken as design variables . 

A method of successive approximations is used to find 
the optimal conditions. 

In the present work, the mathematical progra- 
mming techniques have been used to find the optimal 
design parameters for a single pass rolling mill and 
tandem mill for cold filing. 




1.3 


DESCRIPTIOI'I OF PRESBIOT WORIC 


1*3.1 Finite Eleaent Analysis 

In the present ’'jorli a nathenatical rjcdel is 
proposed for the materiel, rolled in cold rolling 
process using the finite element n-ethod. The mate- 
rial, which is to he deformed during cold rolling, is 
subdivided into triangular elements . Initially the 
material is assumed to be elastic but when it passes 
through the rolls it becomes plastic because of defor- 
mation. Thus the problem is treated as an elastic- 
plastic problem and is modelled accordingly. Depend- 
ing on whether the element is in elastic or plastic 
state, the element stiffness matrix is calculated 
using the elastic or plastic stress -strain matrix and 
the linear strain - displacement relations. The total 
stiffness matrix is then assembled accordiiig to 
standard procedures of structural analysis. By 
applying the boundary conditions for forces and dis- 
placements, one is able to calculate the unknown 
forces and displacements . The stresses and the strains 
in each element can be calculated once we know the 
displacement vector of the continuum. This informa- 
tion about the stresses in each element is used for 
ascertaining whether the element is in elastic or 
plastic state and used for subsequent calculations. 



As incremental displacement method is used, one can 
trace the stiess and the strain history of each element 
as the material passes through the rolls. The informa- 
tion is us-eful to the producer and the consumer of the 
rolled material as mentioned in section 1 . 1 . 

In tiiis analysis one can study the flow of 
material -under the rolls and the extent of deformation 
of the material. Also one gets pressure diagrams for 
the rolls, the general shape of these diagrams is found 
to be in good agreement with the experimental pressure 
diagrams available C25J . 

1.3.2 Design Optimisation 

Further a design problem has been formulated 
as a standard optimisation problem and non-linear 
programming techniques have been applied to solve the 
problem. Two problems have been considered, one in 
which the material is given requisite reduction in a 
single pass rolling mill and the other in which the 
reduction is given in steps in a tandem mill. In the 
first problem two objectives, namely, the minimisation 
of poi^er required for rolling and the minimisation of 
the deformation of rolls, are considered. Constraints 
are placed on the design parameters, which are there 
due to various considerations, like the stresses 




developed in the rolls, location of the neutral point, 
the angle of hite and the deflection of rolls or the 
required power. An upper limit is placed on the 
deflection of rolls when povrer required for rolling is 
taken as the objective function and vice-versa. In the 
second problem, the siim of the defoimiations of the rolls 
of all the rolling mills in tandem, is t alien as the 
objective and the problem is solved such that the solu- 
tion satisfies a prescribed set of constraints . In 
this case, the constraints stated in the case of a 
single pass rolling mill need to be satisfied along with 
an upper bound on the reduction ratio at each mill of 
the tandem. Although the finite element analysis is 
not used to calculate roll forces, torques etc. in 
the present design problem, it can be done ^'/ithout 
much difficulty to predict the behaviour of the rolls 
more accurately. 


*It requires more computer time and storage 




GliAFIER 2 


FIIJITE ELEJIEHT AITALl’SIS OE COLD ROLLING PROCESS 


2,1 INTRODUCTION 

The first step in the finite element 
analysis involves the discretization of the continunm. 
In tiiis study, the prohlem of cold rolling is treated 
as a plane -strain problem. This is justified for 
strip or sheet rolling xfnich involves negligible 
^spread compared to the deformation in the other two 
directions. For a two dimensional problem, one has 
a choice betvreen triangular and quadrilateral ele- 
ments for dividing the continuum. 

The triangular element ia the most obvious 
choice in most of the studies as it is simple and 
versatile for representing arbitrary geometries. 
Further even if one approximates the region by 
quadrilateral or higher order polygons, they can 
always be subdivided i nto a finite number of trian^es , 
Besides, the approach using triangular, rather than 
quadrilateral, elements is much simpler to implement 
computationally [l5j * In a triangular element, 

the displacement function can be assumed to be linear 
as i 


f (x, y) 


a + bz + cy 


( 2 . 1 ) 



wJiereas for a rectangular element, the displacement 
fitnction has to be bilinear as 

f (x, y) = a + bx + cy dxy (2,2) 

Thus all the subsequent analysis and compu*- 
tation is altered because of xy - term. Hence a 
plane constant's train triangular eler'ent is chosen to 
disbretize the continuum. 


2;. 2 1/lRTUAL ¥ORIi FOEinLATIOlI 

The ‘prescriptions' given by Zienkievd.cz 
[l6j for deriving the characteristics of a finite 
element of a continuum is given here in mathematical 


form'. 


■A typical finite element is defined by 


T^odes i.v, 12 and straight line boundaries. Let 

the dris-placemeht of any point vfithin the element be 


’'defined as column vector 




I 


— \ . i'J. • 





h ) 


]h 

f m V 

1 


L V 


(-2^3) 


-in whibh the components of [-Nj are in general 
■ i-Unction of 'position and represents a listing 

’Of -nodal diSpLaoements for a particular element. 



With the displacements knoxm at all points 
within the element, the strains at any point can he 
determined. These will alwa^j^s resiilt in a relation- 
ship which can he written in matrir notation as : 

Assuming an elastic behavior, the relation- 
sliip hetvfeen stresses and strains mil be linear and 
of the form 

[^1 = - {^ol > ( 2 - 5 ) 

6 

Where [ is an elasticity matrix con- 
taining the appropriate material properties ,|6| 
is the vector of actual strains and is the 

vector of initial strains , 

e 

Let define the nodal forces which 

are equivalent statically to the boundary stresses 
and distributed loads on the element ; 



( 2 . 6 ) 



Each of the forces ^ Eiust contain the 

sane number of components as the corresponding nodal 
displacements I and he ordered in -the appropriate, 
corresponding directions. 


The distributed loads 



are defined as 


those acting on a unit volume of material -within the 


element i-jith directions corresponding to those of the 
displacement { f| at that point. 


To make the nodal forces statically equi- 
valent to the actuaJ. boundary stresses and distributed 
loads, the principle of -virtual work is used. Accord- 
ing to this principle, the external and the internal 
works done by the various forces and stresses are 
equated during an arbitrarily imposed (-virtual) dis- 
placements of the nodes . 


[ 6*1 


Let such a "virtual displacements be 
at the nodes. This results by Eq. (2,3) and 
Eq. (2.4-) in forces and strains vjithin the element 
equal to 

and =[b3[®*}' 


respectively. 

The vrork done by the nodal forces is gi-ven 


by 


fjfci® f "V® 

C{8*} ) .(f) 


( 2 . 8 ) 



Siniiarly the internaJL ^/ork per unit Tolume 


done hy the stresses and the distributed forces is 


? ^ f 7 


T 


— 

l* } i 


or ([s*]. - Chi'. { p} ) 


(2.9) 

( 2 . 10 ) 


By equating the external i.'-ork to the total 


internal work, obtained by integrating over the volume 
of the element, one obtains 

( [s*j = ( [s*j " / ( /[b3^ I crjd (vol) 

“ ^ (vol)) (2.11) 

This relation is valid for any value of the 
virtual displacement and hence the equalit3^ of the 
multipliers must exist. On substitution from 
Eq, (2.4) and Eq. (2.5), this leads to 

( y } % ( J [ B]" [oj [_B] d (voD) ^ e]" 

-/[sf - Jw^{ p} 

( 2 . 12 ) 

This elemental relationship can be expressed as : 

[^]= W * Up * '^•'3) 

where p represents the nodal force required to 
balance any distributed loads acting on the element 
and the nodal force required to balance any 



initial strains, such as may be caused by temperature 
changes, if the nodes are not subjected to any dis- 
placements. The first term in the right hand side of 
Eq, (2,13) represents the force induced by the dis- 
placement of the nodes. It can be seen that the 
stiffness matrix is given by 

= j LbI d (vol), (2.14) 

the nodal forces due to the distributed loads by 

d (vol), ( 2 . 15 ) 

and the nodal forces due to the initial strains by 

= -j 1®!^ L E®] a (vol) (2.16) 

The same results can also be arrived at by 
use of the principle of minimum potential energy , ’ 

2.3 ELASTIC PLASTIC AliAlYSlS AS IJOY-LIIIEAK PROBLEM 

In finite element analysis, non-linearities 
occur in two different forms. The first is due to 
material or physical non-linearity and the second is 
due to geometric non-linearity. 

The first category encompasses the p3x>blems 
in which the stresses are not linearly proportional to 



the strains, hut in v/hich only small displacements 
and small strains are considered. I':an 3 ' significant 
engineering problems fall under this ea.tegory and the 
elastic plastic analysis is one of them. In elastic- 
plastic analysis, the elastic stiffness matrix is 
usually repla.ced by plastic stiffness matrix vhen the 
element is transformed from elastic to plastic state, 

2,3,1 Derivation of Elastic Stiffness Matrix 


In this section the stiffness matrix for 
plane const ant -strain triangular element in elastic 


state is derived. 

Fig, (2) shows a 
typical triangular element 
uith nodes i, j, m numbered 
in an anti-clockvdse order. 




are listed as : 




The displacement within an element has to 
he -uniquel^f defined by these six values-. The simplest 
representation is clearly given by two linear poly- 
nomials ast 

u =<i ^ + ^2 ^ - 

. . ( 2 . 19 ) 

The six cohs taiits o^'^, -i = 1, 6, can be 
evaluated easily by salving the t\;o sets of three 
simultaneous equations, which x-ri-ll arise if the nodal 
co-ordinates are inserted '-and the dis place nents equa- 
ted to the appropriate nodal displacements-. One 
finally gets 

u = +tj_-x-+ 0^ y) 

-+ X + C. y) u, + y)u^ 

-V = — ~ (a. + b. X'+’C. y)’v. 

11 11 

(a. + b. -x + C.-y) V. + -(a^ t b^^ X-h y)u^ 


( 2 . 20 ) 



whe re 




"b. = y . - y 


‘=1 = - =^3 


( 2 . 21 ) 


With the other coefficients obtained by a 
cyclic permutation of subscripts in the order i, 3> m 
and 

( 

1 
1 
1 


2 tx = dei. 


^i 

^1 




m 

1 at 

any ; 

its 

thre 


= 2 (Area of 

triangle ijm) 


( 2 . 22 ) 


contribute to the internal uork as : 



3 

J1 

a 

X 

a 

JT 

a 

y 

a 

JX 

"s 

J 


( 2 . 23 ) 


B V 

■8 X 


Using Eq. (2,20) one gets 

r 


1 

2Z1. 


\ 

0 

^3 

0 

m 

0 

0 

°1 

0 

^3 

0 



h 

=3 

^3 






( 2 . 2 ^) 


which defines the matrix £, of Eq.(2.4) explicitly. 



In plain strain, a normal stress cn erists 
in addition to tlie three other stress components - 
For isotropic materials, one 

= E * 

& y = i - iy<- + 

(1. = 2 (1 +iJ) j: 

xy . E xy 

*^2; = i ^ ^ - 

On eliminating a~ and solving for three 
remaining stresses j one gets the matrix as : 

(r^el _ E 

^ J (1 + 1/)(1 - 2 i /) 


The stiffness matrix of the element ijm is 
defined from the general relationship of Eq. (2.14) as: 

01* = jp] W * '"3^ ^2.27) 

where t is the thickness of the element and 
the integration is taken over the area of the triangle. 

If thickrESS t is constant over the element, and, as 
neither of the matrices contains x or y, one gets 

[Kl = t (2.28) 


1-2; I/O 
iJ 1 - 2/0 

1 - 2y 


0 


0 


( 2 . 26 ) 


gets 






>] ) 


a— 

y 


)1- 


(2.25) 


This foim is now sufficiently explicit for computation. 



2.3»2 Derivation cf Fla.stic Stiffness I^trix 


In appl;^rlng the increnental procedure , the 
elastic plastic constitutive lav/ can be expressed in 
terms of finite increments as : 


(2.29) 

The matrix [^D^J is updated for each incre- 
ment by modifying its components because of the change 
in the value of the stresses. Once [^D^J is obtained, 
the element stiffness can be computed by v/riting 

{jf = d (vol) (2.30) 

For a triangular element it becomes 

H ® = LbI ^ Lb^3[b']. t . A (2.31) 


The present formulation relies on the plastic 
stress -strain matrix £dp] for the yielded elements, 
v/here takes the roles of [_D® J . An explicit 

expression for the matrix Von-llises material 

is presented in this section 



The Prandtl-Reuss equation for the devia- 
toric strain increment d ^ ^ during continued load- 

A J 

ing is 9 1 

^^±2 ^ 2G ^ 


(2.32) 



wiie re 


— P ’ _ 

d. X = _ 3 P- ( o '3 "> 

^ 2 * - - 2 * - -/ 

<T^ cr 

Adopting the usual surraation convention the 
equivslent stress cr* and th'e plastic strain increment 

p 

d ^ Sire expressed as folio vs 

^ ( I or^ . ovM , de^ =v^( I dC“~ . d£3/ ) 

-J -J p 

_ (2.^) 

In Eq. (2.33) K = — — corresnonds to the slope of the 

' / _P 

equivalent stresss ( cr) - plastic strain ( J d^ ) 

curve . 


The Von -Uses ;rield criterion and its diffe 
rential form are respectively given by 





(2.35) 


and 

<3r/ d or-/ = fo~do-=^a^ H'^dX (2.36) 

2-0 3 9 

Eliminating d or ^ fro-m Eq. (2.32) and Eo. (2.36), 

2 

2G orj' Cd^''^^ - ^ = 9 

Solving for d and making use of the relation (2.35)? 
one gets 


d> 



/ 



S 




(2.37) 



■with 


a =1 (1 


(2.38) 


Note the equalitj aT:f = erf/ d.^. . in the 

Eq. (2.37), since or/= cr^+ cr''+ o^^is identically 


zero. 


Nov.r substituting d )«. of Eq. (2.37) hack 
into Eq. (2.32) and recalling the definition for the 
deviatoric strain-increment d^r"^. . : 

i J 

/ • - 

d6 . . = dd-.^. - S. . — , d<r.. = ddr + dC- + dC-, 

10 10 10 3 ’ 11 X y z’ 


the deviatoric stress increment d ov. can he expressed as 

3- J 

. q; = 2G 


(2.39) 


The identitjr or^ d^. . = or/ d^lil has been used in 
10 10 


10 10 

the ah ove exnression. 


The total stress increment d of • is, by 

J- J 


definition. 


d or^. - d CJT;. + 3 - 2i/) ®io 


- a /T-'' ^ 2 (1 +^) . G „ . ^ 

~ ^0 3 (1 - 2 2/) ^ii 






SyiiffiBtric 


For a plane-strain problem, it reduces to 


r P', 

[D J = 


E 


t + 


1 -i> 

1 - 2y 

i/ 

1 - 2)J 


Cr— 

X 


/a 


/ I 
CT~ O" 

XY 


- cr 77 

X ^xy 


and d or is given by 

Zjt 


i_^ 5' 

1 -2 2^ " S 
■ 

— o“ 7^ 

y ^xy 


1 

2 


r 


xy 


(2.43) 


^ (1+177 ^ f ”-2> 


C5 CJ^ 
X 2: 


}/ 


) d{lx “^ ( 1-22/ 


+ ( ^ S-.E?.y ) ric- 

+ (. - s ^ ^^xy 


^ / 
cr cr 

^)dC- 

(2.41f) 


Eq, (2.4 i) corresponds to the inverse of the complete 
stress-strain relation l9l of Prandtl-Reuss , The 
original stress-strain relations have been reduced to 
the .single Iq. of (2.4 i) and it must be emphasize that 
the elastic compressibility as veil as strain-hardening 
characteristics of the material are incorporated in 
the matrix • 3y comparing Eq. (2.43) with Eq. 

(2.26), it can be seen that the diagonal elements of 
I^DPjare definitel 3 ?- less than the corresponding diago- 
nal elements of [_D® J. It amounts to an apparent 
decrease of stiffness or rigidity due to plastic 
yielding. 



2,3.3 Criterion for Yielding^ 


A law defining the limit of elasticity under 
any possible combination of stresses is knovrn as the 
criterion for yielding or yield criterion. The two 
simplest and most v/idely used criteria, which are 
consistent with experimental evidence, are the cri- 
terion of Tresca and criterion due to Von-JIises . For 
most metals, Von-1-Q.ses law fits the experimental re- 
sults more closelj!- than Tresca’ s, but Tresca’ s 
criterion is simpler to use in theoretical applica- 
tions, In this work, the Von-ilises yield criterion 
has been utilized. For a general state of stress, 
in terms of the cartesian coordinates, this criterion 
postulates that yielding occurs \^hen 

_ 2 . r 2 2 2 

= 2 LK - "9^ + + (“I - ‘WM 

= 3k2 (2. 

where, k is a parameter dependent on the amount of 
prestrain, k is talten (Y / s/3 ) when Y is the 
yield stress of the material in uniaxial tension. 



2.3.^ Elastic Plastic Behavior of Ilaterials 


The problem is solved for tv/o non-linear 
models of the material. Firstly the problem is solved 
for the elas t-perfectly plastic behavior of the 
material for idiich 


h' = ^ = 0 


d<i 


c p 


(2.46) 


The second model considered in this worli is 
of the tjnpe discussed by itamberg and Osgood £.19 J 
for explaining the tensile stress-strain relations 
for certain materials, Bamberg and Osgood have sho''.vn 
that most of the stress-strain curves in tension can 
be represented by the relation 

e = f (1 + K C ^ ^ (2.4?) 

whereodi 0,02 corresponds to the usual engineering 
definition of yielding. Tensile curves for 
oC = 0,02 are shown in Fig, (4b ) for different values 
of- n, 

H^can be derived from Sq. (2,47) as: 


H 


/ 



n-1 


j; 

cr 


(2.48) 



2.h COIIHJTATIONAL PROCEDUEE 


VJe start mth the strip of naterial divided 
into triangular finite elements as shown in Fig. (3)« 

The stri'p is rolled through the rolls as it passes. 

The boundary conditions are taken as follows : 

F^ = 0, F = 0, external forces on all the 

^ J 

nodes except those which are under the rolls and those 
lying on the line CD, 

F^ = 0, V = 0 on the nodes lying on the j 
line CD, the displacement of the nodes wMch are in 
contact of the rolls are 

u = 0 

V = vertical displacement of nodes because of 
, the displacement AX in X direction 

= A y = -| + H - y^ (2-.if9) 

where 

B^E^= CD +^/R (h^ - h^) - (h^ - h^)^ / ^ +Ax) 

(Xij y^) are coordinates of the node i in 
contact with the rolls and CD as sho\m in Fig, (3) 
is the length of the strip considered for the analysis. 
See Appendix (A) for the derivation of v. 

The computation starts fix)m the stress -free 
condition and the displacement is given to the material 



in X direction in an incremental way. The following 
assiomptions are made in this procedure: 

(i) The arc of contact does not change, i.e,, 
the rolls are rigid. 

(ii) The sprea.d is zero, henco the problem can 
he treated as plane-stra.in problem. 

(iii) To plot the pressure diagram, the radial 
pressure at any point along the arc of 
contact is equal to its vertical component. 

Talcing the above assumptions, the initial non-stead3r 
state of cold rolling is analyzed. 

One thing needs to be noted in this problem 
is that for each node, either the force or the dis- 
placenBnt component is known and hence one solves 
the siMoltaneous equations of the system for the 
unkno>m force and displacement components . 

A flow chart for the computer program deve- 
loped for this analysis is given in Fig. ( 5 ). The 
total stiffness matrix [_k3 is stored in the banded 
form and the nodes are numbered in such a way as to 
obtain the least band width ’B' . A small band width 
is usually obtained simply by numbering nodes along 
the shortest dimension of the continuum, Tlie banded 
form storage of the £k]| matrix requires only II5 



‘ p 

locations as conpared to about N /2 locations for 
storage of the upper or lower triangle of a full 
symnetric matriz of order ’II', 

Also of importance is the savings in 
execution tine that band operations permit. A band- 
form Gauss elimination solver, discussed below^ is 
used to solve the simultaneous equations that result 
after the incorporation of the boundary conditions. 
By -counting the basic operations, i.e., one multi- 
plication or division plus one addition or subtrac- 
tion, the solution time can be estimated to be 
approximately proportional to I'lB /2 in the case of 
band matrix as compared to a value of N'^/b for a 
Xull symmetric matrix. 

The Gauss elimination metliod is applied to solve the 
simultaneous equations. 
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where "u^ and vectors of displacement and force 

components a.re knowi respectively. Firstly we subs- 
titute the value of Icnown in the equations 

corresponding to knoim F2. Thus we get (N - n) 



These simultaneotis 


equations in (l-I - n) unlcnowns U2. 
equations are solved 133^ Gauss band, form elimi ration 


method, the details of i.ixLch are given in CooIc[^17j. 
Once vje Imow the unitnovn vector of modal forces 

F. is calculated from the disnlacement vector 

-V ' f ^1 ? ■ ' 

TJ- '= 1 0:^ r "by substituting the values in first n 


simultaneous equations of Eq. (2.5o), 


Once the nodal displacement vector. is Imo^m, 
one gets the new coordinates of nodaH. points, The^ 
strains and the stresses in the elements can be 
obtained from 

(2. 51) 

and .r i ^ (2,52) 

The listing of the computer program • deve- 
loped for the analysis of cold rolling process. is 
given in Appendix 
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2.5 IlUlEPJGilL EESULTS 
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Hie nuTiiericai resTJlts obtained in the analysis 
of cold rolling process are presented here. The mate- 
rial properties taken for this analysis are as follcn^rs : 

Modnliis of elasticity for the material rolled = 

1,27 X 10^ Tons/sq in. 

Yield strength of the raateidal in tension = ^7*8 Tons/s q in. 
Poisson's ratio = -0»3 

The length of the strip analysed is taken as 
0.8", It is found that strains and the stresses beyond 
this length are negligible, when the strip is rolled. 

The convergence of the results id.th different number of 
finite elements is shown in Appendix C. Two initial 
thicknesses of the strip, namely, 0.064-" and 0.063" are 
considered for the analysis. The final tliickness in 
botli the cases is assumed to be 0,0378.". The roll radius 
is taken as 5" in both the cases , 

Totally four analyses have been conducted by 
using ^00 elements vhth 255 nodes, Tliis corresponds to 
510 degrees of freedom wdth a band-vd-dth of 12 or 
depending on the orientation of the elements , In the 
first two cases, the material is talcen as elastic- 
perfectly plastic with an initial thickness of 0.064”, 



with different orientations of the elerjents as shown 
in Table 2. In the next two cases, the orientation of 
the elements and the initial th-ichness (0.063^') are saiiie, 
but the beiiavioral models of the material were different 
(elastic-perfectly plastic model for the first one and 
the Hamberg Osgood model for the second one). Table 2 
gives the roll forces and torques predicted by the 
present analysis along with experimental values for 
the initial tMckness 0.063” and final thickness 0.0378” 
&i3 ■ Figures 7 and 8 show a comparison of the pressure 
distributions along the roll face obtained for different 
orientations of the elements and for different models 
for the elastic-plastic behaviour of the materials, 
respectively. The figures indicate that the peak values 
of the pressure on the rolls goes up with an increase in 
the roll force. Figure 6 gives the deformation pattern 
and flow of material under the rolls and the locations 
of elastic and plastic zones for the strip having 0.063” 
initial thickness and for a material behaving according 
to Ramberg Osgood, model. 

Table 1 gives tjrpical results obtained by the 
finite element analysis for a material having Hamberg 
Osgood elastic-plastic beha^/iour for a strip of 0.063" 
initial thickness . The columns 2 to 5 of this table 



show the typical results obtained by the computer program, 
whose listing is given in Appendix D. This information 
is used to calculate the total roll force and torque and 
the pressure distilbution along the work-roll interface. 

It can be seen from Table 2 that the numerical 
veiues of the roll force and torque predicted b;'' the 
present analysis are higher than the experimental values. 
The reason for higher values of forces and torques might 
be due to the "rigid rolls" assumption made in this 
work. The elasticity of the rolls can be considered to 
improve the accuracy of the numerical res'uLts using 
Hitchcock's equation, ( 3 . 20 ), Hotjever, this involves 
the use of present anal 3 ''sis program in an iterative loopi 
Thus the total analysis time idii be raiiltiplied by as 
many times as the number of iterations required to achieve 
the desired convergence for the radius of rolls . In the 
present anal^/sis (with rigid roll assumption), the 
execution time was about 3^ minutes on an IBM 7044 
computer. Hence, even if we assume that it takes only 
five or six iterations to attain the desired convergence, 
the total time would be about 3 hours . 



TiSES 1 


Gomponents of Itorce acting at various nodes alcng the 
roll-work interface. 

Initial Thickness of strip = 0.063" Culpit Hiickness of strip = 0*0378" 
Elastic behaviour of liie material ^ Rambeig-Osgood model 


4 25 

3 i 

4 \ 

5 

1 6 

J 7 

iNode No.i 
iSee Sh-g.'l 
1 (6) \ 

F„ 1 

^i 1 
Tons i 

5' 1 

Tons i 

x. 

X ' 

inchi 

I ^i 

1 

i indi 

J Roll Pressui^ 
i at Node i 
i Tons/sq in. 

141 

1.2763 

-1.4670 

0,806136 

0.031067 

-80*5 

146 

-0.0776 

-1.0555 

0.825055 

0.029782 

-56.8 

151 

-0.0012 

-0.9857 

0.843276 

0.0286 10 

-53.8 

156 

0.0484 

-0.9979 

0.861661 

0.027497 

-55.0. 

l6l 

0i281 

-1.0900 

0.880567 

0.0^423 

-59.0 

166 

0.0253 

-1.1751 

0.898623 

0.025464 

-63.7 

171 

0.0193 

-1.2747 

0.917473 

0.024533 

-67.5 

176 

-0.00064 

-1.2935 

0.936393 

0.023670 

-70.0 

181 

-0.0312 

-1.3320 

0.954428 

0.022914 

-72.3 

186 

-0.0514 

-1.3377 

0.973278 

0.022194 

-71.8 

191 

-0.0835 

-1.2944 

0.991660 

0.021560 

-70.5 

196 

-0.0992 

-1.2004 

1.010019 

0.020995 

-66.7 

201 

-0.1217 

-1.1051 

1.027608 

0.020517 

-61.6 

206 

-0.1328 

-0.9809 

1.045805 

0.020087 

-54.8 

211 

-0.1485 

—0 .8268 

1.063535 

0.019736 

-46 .8 

216 

-0.1555 

-0,6592 

1.081144 

0,019442 

-37.4 

221 

-0,1575 

-0.4699 

1.098578 

0.019215 

-27.0 

226 

-0.1518 

-0.2845 

1.115944 

0.019051 

-16.5 

231 

-0.1229 

-0.1423 

1.133030 

0.018947 

- 8.35 

236 

-0.062 

+0.0690 

1.150053 

0.018902 

+ 3.50 


* 

il roll force 


* 

il roll torcjie 



S' = 18.973 -tons 
^i 

S’ (OQ - y.) + '^5' (CQ - X.) = 3.978 Ton- in. 

Xf X ■'x 


e Fig. 3 for notation 




TABLE 2 


Boll forces and torque required to deform iiie material 


Elastic-plastic ji Initial j[E3nal JEI orient j[ Boll X^^orque 

model I Oaiidoiess |Thickness ^configi- | Eorce | Ton-in 

I ! ^ration j Tons | 


1. 

Elastic Perfectly 
Plastic Model 

0,064" 

0,0378" 

1 

I 

5 

18.664 

3«966 

2 


*7 









2. 

- do - 

0.064" 

0.0378" 

1 


6 

17.077 

2.8S2 





I 

! 

2 ^ 


7 



3* 

- do - 

0.063" 

0.0378" 

- do 

16.782 

3.642 

4. 

Bamberg Osgood 
model 

0.063" 

0.0378" 


- do 

- 

ie.S73 

3.973 

5. 

Experimental reailts 





16.56 

2.37 


5 




CHAPTER 3 


OPTIJUSATIOiT OF THE COLD ROLLING PROCE-SS 


3.1 IHTROroCTIOH 

The problem of rolling mill design can be 
formulated as an optimis3.tion problem. Optimisation 
means either minimisation or maximisation of a 
function of the design variables and parameters of 
the system. If there is any restriction on the choice 
of the design variables, the problem is called a 
constrained optimisation problem, othervjise it is 
called an unconstrained optimisation problem. The 
function, with respect to which the optimisation is 
carried out is Imoim as the objective function and 
the restrictions on the design variables due to geo- 
metric necessities are called geometric constraints, 

A general optimisation problem can be 
stated in standard form as follo^-js : 

Find ."x^ wliich minimises f (1^ ) subject 
to the constraints . 

g^. ( x^ ) ^ 0 j = 1, 2, ... m (3.1) 

where is called the design vector and consists 
of the design variables x^ , X 2 , ... f (V-) is 
called the objective function and ( is called 

the j^^ constraint function. 



vsliile designing a rolling mill, the designer 
is confronted with the problem of choosing an anpro- 
priate objective function. The objective can be talien 
as the minimisation of the deflection of rolls or 
power required to drive the -mill. The minimisation 
of the deflection of rolls has to be taken if one 
wants better qualit;r product, while the minimisation 
of pov/er is to be chosen for bringing down the cost 
of p r o due ti on k 

Since both the objectives can not be achie- 
ved simultaneousl3r, one is used as a constraint and. 
the other is treated as the objective in the present 
work. Thus when deflection is. taken as the objective 
the minimisation is carried after deciding about the 
motors available and hence their output vovjev puts 
an upper bound on the power required for the rolling 
process. In the second case, the minimisation of the 
powder required for rolling is carried by placing an 
upper bound on the maximum possible deflection of the 
rolls , 

In order to evaluate the objective func- 
tions and the constraints , the expressions for the 
roll force and torque are required in terns of the 
roll dimensions , the exit and the entry tensions and 
the dimensions of the material being rolled. These 
expressions are derived in the following sections. 



3.2 EXPEESSIOIIS FOR EiE OBJECTIVE AID 
CONSTPAIWT FUIICTIONS 

In practicej for cold, rolling, the strip 
has to xoass through three zones. 

In the first zone, i.e., 
near the plane of entry to 
the roll gap, the strip is 
elastically compressed 
until the yield stress is 
reached, in the second zone 
it is further compressed 
plastically and in the 
third zone it recovers 
elastically, Fig. (9 ). 

3.2.1 Elastic Deformation of a Strip 

By considering the force balance of a 
transverse slice of material rolled, shorn in 
Fig. (1), one gets the basic differential equation. 

= 2-^ R (sin 0 + /Ul cos 0) (3 •2) 

Bland and Ford made the following assump^ 
tions to derive the normal roll pleasure from this 
basic differential equation for the rolling process. 


Zone 1 



ELG. (9) 



(l) The a,rc of contact remains circular even when 
roll flattening occurs . 

constant for the v/hole length of the arc 
of contact , 

(3) The elastic compression is negligible (This 
assumption is waved late r_ when the elastic 
deformation is also taken into account), 

(4-) The spread is zero, 

(5) Von-!^Iises criterion of plasticity is employed, 

(6) The deformation is wholly- homogeneous , 

(7) The radial pressure *s’ at any point along the 
arc of contact is equal to its vertical component. 


The normal roll pressure is given by [63 : 


r h crj- /UH 


(3.3) 


from exit to the neutral point and 


, , OT - H) 

from entrj'" to the neutral point, where 




H = 2’^^ tan"”^ 0) (3.5) 

and 

= 2 '^^ tan-l C ^ ) '3.6 ) 



(3.7) 


0 ^ is given ■b:^ 


0 . 


' aur ti t i m 

J (li^ •’ II2) 

iTT 


The Toll force ?_s ooiainod. from "tlie expression 

/ 2 " 

= R / d 0 + R* y s“ d 0 

o 


(3.8) 


0 . 


n 


where 0^ ,, the Single of Lhe neutral point is given by 



The roll boreue is given by the expression 


G 



s 0 d 0 + 


(sTj h, - ajij) j 

J 


(3.11) 


In. these equations j the tensions and denote 

the absolute values , the signs having been taken care 
of in the deri’'mtiono 


The problem in evaluating P from equation , 
( 3 . 8 ) lies in the fact that the yield stress k occurs 
as a variable in equations (3*3) and (3.^)* Ford and 
Bland ^6j has -shorn that vhen mean yield stress is 
used in place of ns’ as an approximation the error mil 
not be more than 20 for cold rolling. 



WhBn the mean 5d.eld stress k is taken 8ts 
a constant, the e:xpressions for the roll force and 
the torque can he given in terms of certain non-^ 
dimensional quantities, as 

/ ' ^ 

S = k E (h^ - hj) . (1 . J ) . f (a, r, b) (3.12) 

k ^ 

and 

G = R k (h.j - h 2 ) . (1 - U) . (a^ r, b) (3^13) 

k 

li^ere non-dimensional parameters are given by 







b = ( 1. 

k k 

and r = ( li.^ - h 2 ) / 


The functions f^ and 
tion of the expressions for P and 
Appendix (B) , 


and the deriva- 
G are given in 


3,2*2 Elastic Deformation of a Strip 


The effect of elastic deformation in cold 
rolling of strip is three fold, the total roll force 
is thereby directly'’ augmented; the friction between 
the roll and the strip over the elastic arc increases 
the "friction hill" there by increasing the plastic 
roll force; and the effective roll radius is maiified. 



Except for -very small passes, tiie elastic 
compression at tlic entry causes negLigilDle difference 
in the roll force. The contribution of the elastic 
effect at exit, however cannot be neglected; in small 
passes it may account to of the total roll force, 
and in normal reductions, it is seldom less than 
[ 21 ] . 

Mailing certain assumptions and approxima- 
tions,. Ford Ellis and Bland [^22_^have shoTO that the 
elastic contribution to the roll force is given by: 



( 3 . 1 ^) 


( 3 . 15 ) 


( 3 . 16 ) 


( 3 . 17 ) 



and bhe effective front tension at the "beginning of 
the plastic zone is 


Oj' = OJ-- 


h. 


( 3 . 18 ) 


niinimum thickness of the strip in the 
pass, it is given "by 

h^ (1 - (k - err) 


h = h h cr* + 


E 


(3.19) 


The Hitchcock’s formula for the radius of the contact 
arc l/, naj'iiely, 


{ 16 (1 -^^) 

R = R (1 + ^ 


Tf E 


• h^ - h2 


( 3 . 20 ) 


is modified, due to the presence of elastic zone, 

/ 16 (1 -' 02 ) :p 

R = R (1 


Tf 


E 


(/6 + + 6 ^; + /Brj 


(3.21) 


where 


6 = h^ - h2 

h^ ( 1 -9^) 

^2 = S 

(1 - 0 ) 


( ka - d^) 


and = 


E 


(ha ** d:f) 


Thus the total roll force and the roll 
torque, hy considering the elastic deformation, 



are given by 



( 3 . 22 ) 


and G = Gp + G^_^ + G^^ (3-23) 

where and G^ are roll force and toroue due to 
plastic deformation of the strip. 


To cal cula';e the foree and the torque due to 
plastic deforma'i on, i-je ; se the new effec’ ive back and 
front tensions given by oquation (3.17) and (3.l8), 
Also we use h^,, given by equation (3.19); instead of 
h 2 SIS the plastic zone spreads upto the r'lnimum thick^ 
ness between the rolls, beyond which the elastic 
recoveiqr starts. 


To caloilats the exact value of the roll 
force and the torque an iterative procedure is used. 
The equations for the roll force and the torque in- 
volve the modified radius R^. When calculation is 

/ 

started R is assumed to be equal to R, Once the roll 
force and the torque are calculated, the calculated 
value of the roll force is used to find R from modi- 
fied Hitchcock's formrla (3.21). The new value of R^ 
is now used to calculate a more accurate roll force 
and torque. The procedure is continued till the 
desired accuracy for the roll force or the roll radius 
or for the both is achieved. 



3.2.3 2617313 511318 Loads for Rolls 


Once the torque required to roll the material 
and the roll force are known, the stresses present In 
the rolls can he calculated^— The upper hoiind on these 
stresses Is going to a constraint In the optimisation 


problem. 

With respect to 
Fig, ( 10), Wusatowsld. [^23j 
gives expressions for the 
bending stresses, under 
the action of roll forces 
and the torsional stresses 
due to roll torque. 



FIG. (10) 


The bending stress is given by 




32 _ 

TfD^ 


b. 




(3.2^) 


where ■= b + C. 

If the material is rolled in the centre of 
rolls, which is the case assumed in the present work, 


b. P b. 

^ = 7fi3 ^ b - - > 


(3.25) 


The torsional stress in the roll barrel is given by 


2 

tfE^ 


T 


( 3 . 26 ) 



The expression for the compound, stress under . 
the simultaneous action of bending and torsion is 
given by 

/ ; 

erg - “V steel rolls C3.27) 

and = .375 + .6^5 2 cast- (3.28) 

iron rolls , 

3 . 2 .4 Deflection of Roll s 


Once the roll force is known the extent of 
deflection at any point froii the end of barrel to 
midway along the length of the roll can be found foTom 
the equations 


n 

.. P (x - 

yB 

- 2^~1' 

and C o 

5 P 

ys 

TfE D^ 


b. 

where n 

= Id •• — ^ 


n 


2 b^ n - 2x (x + 
(3.29) 


L 


(2b ~ x) - n (2b -- n) 


] 


(3.30) 


the meaning of the other s-ymbols is indicated in 
Fig. (10). The deflection at the ends of the roll 
barrel (at x = n, and x = 2b - n) is equal to zero. 
Hence maximal deflection tri-ll occur at x = b = L^/2. 


As suffix B denotes the deflection due to 

S due to shearing forces, the total 


bending forces and 



(3.31) 


deflection of rolls is 'given by' 

" ^yB ^ S’-s 

3.2.5 Energy Co n sumed in Gold Rollins 

Broadly speaking^ the total energy required 
for rolling can be divided into four parts: 

(1) the energy needed to deform the material as 
it passes through the rolls, called basic 
energy of rolling, 

(2) the energy needed to overcome the frictional 
forces in the roll neck bearing, 

(3) the energy required to balance the frictional 
losses in pinions, reduction gear etc., 

(V) the energj/' required to compensate for losses 
in the various electric motors, generators 
and the complete electrical circuit. 

Further when significant amounts of coiler, 
decoiler tensions are employed, some energy would be 
needed" to pull the strip through the rolls. In this 
work, the energy consumed in rolling is computed by 
Gonsidiering parts (1) and (2) alongwith the energy 
due to coiler an'^ deeoile’^ tensions. 



The basic energy of rolls is nothing but the 

■work done in overcoEiing the torque required to roll 
the material, where 


work done per roll = 2~(f G per revolution. 

Hence basic energy for two rolls = kJfG per 
revolution . 


The work required to overcome the frictional losses 

P / 

in the roll neck bearings of two rolls = 2 X 271 ^ /A. d 

= 2 TT P /-V d 

The coiler and decoiler energy consumption per minute 
is given by the following expressions: 


Energy consumed per rnimte (coiler) = F v 
Energy'- consumed per minute (decoiler) = B v 


In 

h. 


where P is the pull and B is the drag and v is the 
speed of strip being coiled on the coiler. 

If b^ denotes the width of the strip, 

F 


Hence, 

Energy consuraed per minute (coiler) = ^2 "^1 * ^ 

Energy consumed per minute (decoiler) — cy!j*’h 2 ^ 



Total energy consumed by coiler and decoiler 


= ( cr^) Ij per revolution 

where, rl is the number of revolutions per minute of 
the rolls . ■ — 

Total energy consumed per revolution for strip rolling 
in the mill 

/ hp b V 

= ifir G + 2l1 P/-«- d + (d:| + oj) 

h^ b V , ^ 

Here ““jl — can be expressed in terms of the other 

dimensions of the roll as 

= sir E r 1 + ^ (1 - cos J (3.32) 

where, 0^ is the neutral angle and is given by 
Bq. (3.9). 

Therefore, total energy consumption in cold rolling 
of strip In the mill 

= ZtjjQ + d + R ((J3 + (1 + ^ C1 - cos 0^))J 

(3.33) 

3.3 COHSTPuOTT SET 

For the optimisation of the objective func- 
tions considered, the factors that impose restrictions 

on the design parameters need to be considered. The 

I i,-. r .'"t;FUR 
i./.L J jB RARl 



restrictions come from variolas considerations like the 
stresses developed in the rolls, the location of the 
neutr^ point, the angle of bite and the deflection of 
rolls / the required power. The following constraints 
are considered in the present formulation: 

(1) The lower bound on the radius of rolls is taken 
as zero, ho upper bound is considered on the 
radius of the rolls as the power restriction 
(either as a constraint or as an objective) 
limits its value. 

(2) The lovrer bound on the back and front tensions 
is taken as zero. The upper bound is not con- 
sidered because of the same reason as given in 
( 1 ). 

(3) If the neutral angle, 0^, is less than zero, 
the rolls will skid over the material and if 
it is greater than the angle of bite, 0^ ; the 
rolls jam and will not move. Hence 0^ is 
restricted as 

0 0n ^ 

(2+) The angle of bite, 0^, is constrained to be 
positive as 

The may imuTn angle of bite is dictated by the 



friction between the rolls and the material* If 
the angle of bite is more than twice the an^e 
of friction, we need to apply front tension. 
Hence, 0^ is limited as 

0^ ^ 2 tan"'* ifUH.) 


Hic compound stress in the rolls due to bonding 
and shearing forces should not be more than a 
certain permissible value, which is dictated by 
the material of the rolls : 



Wlien minimization of the deflection of rolls is 

taken as the objective, the power required for 

the rolling process is restricted as; 

/' 

2f|( 2G + + <5^)(1 + ^"5 - cos 0^)) 

where P is the minimum power available, 
max 

If minimization of power is taken as the 
objective, the deflection of rolls is limited to 
a maximum value of 3-s • 

Cy = Cyg + Cyg 4 

where € is the maximum permissible deflec- 
msix 

tion, which will be decided by the quality 
requirements of the rolled product. 



3.^ STAISICNT OF THE OPTIMISATION PROBLEM 
FOR A SINGLE PASS ROLLING 


For a single pass rolling, 4;iie design vector 


is taken as 


X 



W ©2 J 


( 3 .^) 


Tlie expressions for the objective fwiction 
is given by ' . 

f . ( t ) = c.^ = Cyj + Cyg (3.35) 

•where and C are given by eq-uations (3 •29) 

and (3,30) respectively, if deflection is taken as 
the objective,. 


or 



if powei Is taken as the objective. 


The nine constraints can be stated as t 

(1) - R ^ 0 

(2) - sf 0 

(3) 

( 4 ) - ^ 0 

where is given by Eq, (3»9) 


(3.36) 


(3.37) 

(3.38) 

(3.39) 
(3.^0) 



( 5 ) 

■0n-^-r<o 

( 3 *^ 1 ) 


where 0^ is given by- Eq, (3,7) 


(6) 

>• 01 ^ 0 

( 3 . 42 ) 

( 7 ) 

01 -^tan"'’ i^) ^ 0 

( 3 i 43 ) 

(8) 

^ ” ‘permissible ^ 

(3.43a) 


where or: is given by Eq. (3.27) 



(9) • 2Tf r 2 G + R (<7r+ (T^) , 

W f * 

(1 -cos V)]- ° (3.>4.) 


when deflection is taken as objective, 






- e 


i ° 


(3A?) 


where and Cy.g are given by equations 

(3«29), (3 ‘SO) respectively, when power is 


talcen as. objective. 


3,5 STATll’IEUT OF TiiE OPTIMISATION PROBLEM POE 
lULTIP-^SS ROLLING IN TA 1 ®EI 4 KELL 

In the case of material being rolled in a 
tandem mill ha-^/ing m single pass rolling mills in • 
series, the objective function is taken as the sum of 
deflections of t]\e rolls in all the ’m’ rolling mills. 



Hence the ohjective function in tandem mm 
rolling is given by: 


i=1 ^ ^ 

where ' and are the deflections due to 
bending and shear stresses in the i"*^^ mill and are 
given by equations ( 3 . 29 ), ( 3 . 30 ) respectively. 


i3M) 


The design vector, 



is taken as 



(3.^7) 


't/ll, 

where is the radius of the rolls in the i mill, 

cnj is the back tension in the first mill, is the 

front tension of the last (m ) mill, is the ten- 
sion between i - and i"*^^ mills (i = 2, 3, ... m), 

ttl 

is the thickness of the strip between i - 1 and 
i^^ mills (i = 2, 3, m), is the input thictoess 

and is the output thickness of the strip. 



The constraints can he expressed as : 


( 1 ) 

- h i ° 

i — ”1 y 0 . . m 

(3.if8) 

( 2 ) 

4 0 

i = 1 , m + 1 

(3.if9) 

(3) 

-^nl 4 ° 

i = 1 , ... m 

( 3 . 50 ) 


where" 0 ^ 

til 

is the neutral angle in the i 



mill and is 

given hy Eq, (3-9) • 


(if) 

«>nl - 

0 i“* 1 ?.**m 

(3-51) 


where 0^^ 

th • 

is the an^e of hite in the i 



mill and is 

given hy Eq. (3.7). 


(?) 

“ ^ 1 i ° 

i = 1 , . m 

( 3 . 52 ) 

(6) 

0^^ - 2 tan 

(/*o^)s^ 0 i = 1, ... m 

(3.53^ 


where^it*-^ is 

the friction between material and 


■feh. 

rolls in the' i mill. 



(•7) 




cur. ' -* 

■ permissihle 

where <S^ 

•^permissihle 


the total 
. th 


( 8 ) 


• 0 ■ i = 1, ••• m (3*?^-) 

is maximtsm peimissihle. 
stress for the rolls?, 

stress on the rolls in rolling in the i '' mill 
and is given hy Eq. (3*27). 

2ir[2 Gi + Pi/A-^1 di + % 'X "■ • 

— (1 - 004 ° 

1 = 1 , ... m C3-55) 


(1 + 





3,6 SOLUTION TECHI'IIQUE 


The optimisation problem formiilated in 
sections (3»^) and (3*5) can be seen to be of the form 
given by equation (3,1). Keeping in mind that an 
unconstrained optimisation problem is alv/ays easier to 
solve compared to a constrained problem, the method of 
Sequential Unconstrained l-iinimisation Teclinique (SUMT) 
has been applied to solve the constrained minimisation 
problems formulated earlier. The steps involved in 
finding the solution by SUM! method are as follows : 


Step 1 


Determine an Initial feasible vector , 

which will satisfy all the constraints, findirc^ 

'x^ , care- has to be taken to see that it will be a 
o ’ 

free vector but not a bound vector. A free vector 
is a point which will satisfy all the constraints 
with strict inequality and a bound vector is one 
which will satisfy at least one constraint "\/ith 


equality. The reason for avoiding a bound vector is 
that if one or more constraints aia satisfied with 


strict equality sign, the term 


m 

ZL 

3=1 


h 




which will be needed in Lhe subsequent steps, will 
become infinite. 



Construct a ne\r function 0 as 

ra 

0 i'x ,Y*) = a. ( tp-)-- (3.57) 

3=1 n 

where ’T' is some positive constant. 


In the present work, the starting value of 
is chosen so that the two terms in the expression 
for 0 are approximately equal when evaluated at the 

point . 


Starting value of 


f ( 


^o ^ 


m 

r 

3=1 




(3.58) 


step . 3 

Starting from the vector x^ , find the 
unconstrained minimum of 0 (?) for a fixed value of 
Y^= y/ . The Davidon Fletcher Powell (DFP) method is 

used to minimise the 0 function. 


This method involves the fell- owing steps . 


( 1 ) 


Start with a given vector x^ and an initial 
positive definite symmetric matrix . 

can he taken as a- leintity matrix. 

Then set S^ = - [ ^o j ^ ^o 


(3.59) 



where 0^ is the gradient vector of the fnnc- 

3 .nd S _ is cs-llsd 




it ion ^0 BVcilns-t^d s-t x 
the direction vector. Set the index q - 0 for 
generalizing the procedure. 


(ii) fizid a new vector x^^^ as 


X 


■q+1 


= * ‘^a 


(3.60) 




where «; ^ minimises 0 ) . 

^ q* is kno™ as the optimum step length along 
the direction 7^ and step (k) gives a prooe- 

dure for fiudingix;^ . 

. t is the gradient vector 

iii) Evaluate |7 0^+^ j ’ 

npw search 


, rt ^t X . Generate the new search 
of 0 at x^.j,.| • 

direction as 

si - = -LViJ '^^“-+1 

where = [h,] t[M^] 


(3.70) 


with [_ J ” K q 


s . s 

2- 

^ • h 




r Y, 1 = y ^ti ' ^4 




LM 


and 



(iv) Onc8 tliG HGw SG3,rcli dirGction S . is known 

Q't‘1 ^ 

set q =: q -f 1 and repeat steps (ii) tlirougli 
(iv) untill the function jZ)' is iiinimised. To 
f ino. . v/liether the function j3 is minimised in 
. any: particular iteration q, a convergence 
„ , criterion has to Ice used at the end of step 
(ii),. . ‘ ' 


Step 4 


This step indicates a method for. determining 
the optimum step length v ^ , from a given point and 

along a particular search direction. The cubic inter- 
polation techniques has been applied for this purpose. 

, ' is the quantity uhich minimises 

0 ( + 5 <q^ '"'q where and are already 


knovm. Hence 0 ( x' +p(„ S ) can be treated as 

H M, m 

a function alone and it is known as the one 

dimensional minimisation problem. 


Vlien subscript q is dropped, let 
f (<?<) = 0 ("?•!• .o( "^ ) 


( 3 . 71 ) 


Apr r ox j mating f (c< ) by a cubic polynomial 


h (o<'),. f ( 5<) can be written as 

f (;x^ ) -Izr'h ( ^ ) = a .+ bp( -f- c»<^ + d^^ 


( 3 . 72 ) 



Let for A = 0 
and for CK' = B 


f = f^ and / = f y ^0 

/ J 

f = fg and f = fg 7 0 

Then the value of cs,/ which minimises h (K) is 
given by 


o< 


= A 


f. + 2 -!* Q 


f,/ f-o ZZ 


( B - A) 


(3.73) 


where Z = 


Tb" 


.o N 

l.t3_ 

T" 


/ / 

+ f 'I- f 

A 


/ /Ji 

■P P ) ^ 

"A '"B^ 


(3.74) 


and Q = (Z"~ - 

After findfng eC , convergence criterion 

llh,AL /f 

where V = should be checked to see whether 

also coireivyonds to minimum of f (p<^). If the 
convergence crite:.:ion is satisfied, (;( can be taken 
as '' 5 otherwise the procedure has to be repeated 
as stated below. 


d'" . / * 

Bc-aluatc ar 


t/D 


It will be either 


positive or rn i',:at‘'.ve , 
df 1 


i) If :s pc-itiye, set B = ^ , 

f ( e< '■'} and fg^ = ^ ( ^'^•'and calculate 


"B 


a new value of using equation (3.73) » 



negative, set A = o( 


'A 


f ( «( *) and 

# 


df ! 

^ calculate 


a new value of ci( as per equation (3,73). 

Once a new value of < * is found, test for 
convergence according to equation (3.74) and the 
subsequent procedure has already been stated above. 


The steps 1 to 4 indicate how the function 
0 is minimised for a particular value of y*" . At the 
end of the minimisation process, a new vector, which 
minimises 0 ( 1?, yy) , will be obtained. IIow the com- 
plete procedure is repeated with this new vector as 
the starting point and with a value of smaller 
than the previous value , The process is oontinued 
for several values of and can be terminated when 
some or a.11 the following conditions are satisfied 
simultaneously'. 

(1) Final value of is very small as compared to 

(2) . for two conse qutiveV^is are almost equal, 
■^min 

(3) The value of f corresponding to 0^^^ for two 

consequtive are almost equal, 

(4) There is no significant difference between the 
values of 0. and f. 

The design vector corresponding to 0^^^ for 
two consequtive K' s remains almost identical. 


( 5 ) 



3.7 IJUilEEICAL RESULTS FOR SIRGLE PASS ROLLIHG 


In the present work, oast steel is taken as 
the material of rolls and 0,08^ carbon steel is talien 
as the material to be rolled in the rolling process. 

Numerica2 Data; 

IfeLxirnm permissible stress in the rolls = 637 tons/sq in. 
Modulus of elasticity of the roll material = 1 .275 x 10^ 

tons/sq in. 

Coefficient of friction for roll neck bearing = 0.08 
Constant in Hitchcock's equation (c) = 1,67 x 10*^ 

Modulus of elasticity of the material to be rolled 

k 

= 1.27 X 10 tons/sq in. 

Poisson’s ratio for the material to be rolled = 0.3 
Coefficient of friction beti-Teen the roll and the material 

= 0.055 

Average 3d.eld strength of the material to be rolled 

= 37-6 tons/sq in. 

The actual yield strengths used for different 
reductions are given in Table 3 l21j. 

TABLE 3 

Actual yield strength for different reductions 
(for 0,08^ C steel) 

^■■Heduction | 0 ^ | 10/ | 20^| 30 ^ | 40^ 50 ^ | 6C^ | 70^ 

30.?g 37i^^1.2,^M..o!^6.2g.^8.og ■^9-7 



3.7 *1 Dsil scbion of Rolls as Ob.iectiv© 
f or I'-tlnimisation 

The design optimisation problem is solved for 
three inj^ut thicknesses ( 0 . 13''3 0*063" , 0.013") >ri.th 
four percentage reductions for each thickness (lO^, 20/1, 
30 /^, 4-0^0. Each of these twelve problems is solved by 
using three different upper bounds on the required poi-zer. 
The numerical res-ults are given in Tables 4 to 6. The 
following trends can be observed from these tables, 

(i) The stresses induced in the rolls have been found 
‘ to be ver 3 ^ small, because in rolling of thin 

strips, forces ererted on the rolls are small, 
and are of the order of 17 tons in the case of 
strips of 0.063” Initial thickness having 40^ 
reduction. 

(ii) The constraint on vo\ier has been observed to be 
active in all the problems. 

(iii) For very thin (0,013”) strips, the back and the 
front tensions at the optimum point are consi- 
derably higher compared to those of thicker 
strips (0.063”5 0.013”)* 

(iv) As the percentage reduction increased, the radius 
of rolls slightly decreased and the tensions 
increased in the case of very thin strips. In 
the case of thicker strips, all the three 



design variables are reduced as the jercentage 
reduction is increased. 

(v) In strip rolling the power required to deform the 
material is supplied through the toraue apulied 
across the rolls and the back and ; front tensions. 

In general tjie power available i-ri-ll be distributed 
betx«reen the torque applied tlirough rolls and the 
tensions It has been observed that in the case 
of very tliin strips, an increase in the available 
power (because of a higher value of P ) resiilts 
in increased tensions and the radius of the rolls 
do not change much (rather decreases slightly). 

This means that the work done by tensions is more 
and the torque on the rolls \i±ll be less, which 
also implies that roll forces are reduced result- - 
ing in a smaller deflection of the rolls. 

(vi) ■ In the case of thicker strips, the optimum tensions 

do not change considerably but the roll radius 
increases which implies that more power is required 
for rolling. With an increase in the roll radius, 
the stiffness of the rolls increases resulting in 
a smaller deflection of the rolls . 

(vii) In Figures 11 to 13, the behavior of the objective 
function is shoirn plotted against the percentage 
reduction for different values of the maximum 



power avail abe to the rolling mill. It can be 
observed from these figures that for very thin 
strips, the availability of more power has a 
negligible effect on the objective function. 
However, in the case of thicker strips, the avai- 
lability of more power results in much greater 
reductions in the deflection of rolls. 

3,7.2 Power Required in Rolling: as The 
Objective For Hinimisation 

In this case, the design optimisation problem 
is solved for one input thickness (0.063") with four 
percentage reductions, namely, 10;^, 20 %, 30 % and 40 ^. 

The deflection of rolls is taken as one of the constraints 
for the problem and an upper limit of 5*0 x 10 inch is 
.put on it. The results are given in Table 7 . The trends 
observed in this case are as follows : 

(i) The stresses in this case also are very small, 
because of small roll forces. 

(ii) The lower limits on the back and front tension 
becaiiie the active constraints at the optimum 
point. In fact, the back and fron tensions are 
too small for all practical purposes. 

(iii) It can be seen that a more efficient way of 
reducing the power is to decrease the radius of 
rolls, i.e., to reduce the torque applied 
through the rolls. 



3,8 rlUlMPJG/'iL RESULTS FOR mi TI PASS ROLLI'^G 
III TAHDEI'f MILL 


The material properties for the rolls and the 
material rolled are same as those are given in section 
3.7* The problem is solved for a tandem mill having two 
single pass rolling mills in series. The maximum reduc- 
tions aH.lowable in the first and the second mill are hOfs 
and respectively. The maxirAum allowable powers for 

both the mills are same. It is 800 ton-in per revolu- 
tion for both the nhlLls . The total reduction given in 
both the passes is ^ 0 % and the initial tliiclmess of the 
strip is .13”. 


The initial and the final design vectors are 


given belo'w; 



V. 00 000 


|^)+. 00062 ] 

If .00000 

i 

if .00047 

0.07000 

. (v 1 - ' 

0.07003 ^ 

0.07000 

! ’ 1 1 

: 0.07003 

0.07000 


0.07001 

0.087 

1 

0.078036] 

4 , 


Obisctive function at the start of 

op tiniis ati on - 

Optimal objective function " 


.7253 X 10“^ in. 
.7055 X 10“^ in. 



Tlie total nmler of one -directional optimisa- 
tjon steps are 21 and the total computer time talcen is 
^8 minutes on lEil 70^i-^ computer. The observations are: 

i) the tensions and the roll radii are almost same 
at the optimum point as at the starting point, 

and 

ii) the first mill works to its highest capacity, 

that' is the maximum possible reduction has been 
obtained in the first pass itself. 

The constraint v/hich became active at the optimum point 

vas the upper limit on the maxim™ reduction allouahle 

in the first pass. 
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TABLE ( 4 ) 

Initial Thickness of the Strip Being Boiled = ,13” 


T 

I* 1 

- TTIQ-V 

[ Design Vector 

1 

Obi . fn • fen wer Tifo 

- of iCciraiter 

! ^ • I 

lEeduc- 

|tion 

r Ton- in/' 
|rev. 

1 X. 

1 ^ 

1 

1 Inch,. 

i X 2 t x_ 

f ^ f 5 < 

|Ton/ |Ton/ | 

|Sq. ln.| Sq. ini 

(detlecjiused Jitera- 
-tion) 1 Ton-in/itions 

-4 1 1 

10 in.t f 

Itime 

fsiin :Sec 



Initial 

3.000 

0.384 

0.427 

4.372 





800.00 

Optimal 

8.102 

1.024 

1.104 

2,992 

781.04 

25 

09 I 52 

20^ 


Initial 

3.000 

0.385 . 

0.481 ■ 

9.832 

mm 



800,00 

Optimal 

6.239 

0.622 

0.614 

7.459 

794.518 

11 

05 : 09 


305 S 


Initial 

3.500 

0.385 

0.550 

42.33 




800,00 

Optimal 

5,199 

0.253 

0.278 

13.94 

794.32 

23 

05 t 01 




Initial 

4.000 

0.385- 

0.641 

31.260 




40^ 

800.00 







14 

05 { 01 



Optimal 

4.419 

. 0.463 • 

0,644 

23.660 

795.49 



10^ 


Initial 

8,063 

0.429 

0.293 

3.032 




900.00 

Optimal 

8.711 

0.436 

0.081 

2.565 

854.2 

34 

11.: 




Initial 

6.239 

0.622 

0.614 

7.459 




20$fa 

900.00 





892. 14 

26 

04 ; 55 



Optimal 

6.604 

0.531 

0.402 

6.498 





Initial 

5. -199 

0.253 

0.278 

13.94 

- 



30fci 

900.00 

Optimal 

5.636 

0.370 

0.148 

uHl 

887.-^ 

23 

05 : 11 



Initial 

4.419 

0.463 

0.644 

23.660 

- 



40^ 

900.00 





897.19 

30 

10 : 1‘‘f 



Optinal 

_4.907 

0.233 

0.352 

17.833 





Initial 

8.711 

0.436 

0.081 

2.565 

- 


14 l 00 


1000.00 



2.104 

972.6 

17 


Optimal 

_8.716 

0.631 . 

0.078 


— 

20^ 

1000,00 

Initial 

6,604 

0.531 

0.402 

6.496 

— 

12 

05 i 06 


Optimal 

7.304 

0.555 

0.977 

5.124 

996.32 




30$^ 


Initial 

5.638 

0.370 

0.148 

11.300 

— 

15 

05 : 06 

1000.00 




9.6^ 

977.09 


Optimal 

5.987 

0.853 , 

0.590 


.. 

40% 


Initial 

4.907 

0.233 

0.352 

17.833 

— 

22 

07 : 11 

1000.00 

Optimal 

5.285 

0.368 

0.381 

14.674 

993 .83_ 




TJSBIE (5) 

Initial- Eiickness of the Strip Being Boiled = ,063” 


TTp 


I 

iEednc-JI 
:tion 


max 


frev. 


i Design. Vector 


| .x, j X, j 

[ Inch, j Ton/ ] 
i 1 Scb to 1 

1 

[ Ton/ 

( Sq.in 


Jobj.fh.i Power | No. of| Coai*u- 
][(defleci[ used | Itera-| ter 
I -tioa)| lon-im ticos f IMme 

1 t I Min: Sec, 




Initial 

3.000 

0.794 

0.882 

32.590 

mm 




10^ 

800.00 

Optimal 

9.442 

2.464 

2.900 

2.804 

785.64 

12 

05 : 1 

0.. 



Initial 

3.000 

0.794 

0.992 

42.160 

mm- 




20% 

800.00 

Optimal 

7.723 

2.082 

2.416 

3.580 

769.28 

65 

11 : 

56 

-/ 


Initial 

3.000' 

0.794 

1.133 

51.650 





30% 

800.00 

Optimal 

6.895 

1.055 

1.340 

5.544 

783.472 

52 

12 ; 

1? 



Initial 

3.000 

0.794 

1.322 

57.03 

.. 




40^ 

800.00 

Optimal 

6.117 

0.752 

1.170 

8.280 

786.4 

51 

12 : 

05 



Initial 

9.442 

2.464 

2.900 

2,804 

- 




10% 

900.00 

Optimial 

10.242 

2.320 

1.613 

1.476 

880.84 

21 

05 ! 

03 



Initial 

7.723 

2.082 

2;416 

3.580 





20% 

900.00 

Optimal 

8.371 

2.257 

2.638 

3.019 

874.69 

20 

05 : 

1C 

30% 


Initial 

6,895 

1.055 

1.340 

5.544 

- 

28 

10 ! 

09 

900.00 

Optimal 

7,430 

1.039 

1.714 

4,700 

860.59 






Initial 

6.895 

0,752 

1.170 

8,280 

- 


06 : 

50 

40% 

900.00 

_ Optimal 

0.770 

1.228 

-i-oSlt 

885.^ 

20 

6.632_ 









10% 

1000.00 

Initial 

9.442 

2.464 

2,900 

2.804 

f.* 

41 

12 : 

07 


Optimal 

11.560 

2.373 

2.560 

1-.185 

976.31 





20%. 

1000.00 

Initial 

Optimal 

7.723 

9,156 

2.082 

1.742 

2.416 

2,103 

3*580 

2.524 

961.47 

62 

12 s 

OT 

30% 


Initial 

6,895 

1.055 

1.340 

5.544 

- 

41 

13 : 

56 

1000,00 

Optimal 

8,045 

0-998 

1.378 

3.972 

9^3.64 






.40% 


Initial 

■ 6.117 

0.752 

1.170 

8.^ 

- 

33 

14 i 

07 

1000,00 

Optii^ 

7. 102 

0.827, 

1.382 

5.903 

983.93 


— 

— 



Initial !I!hio]mess of Strip Being Soiled = ,013" 



i P \ ' ■ 

DesiOT- Vector 





jReduc-, 

' lUCAJi. V % 

f Ton- In/v 1 

J T’O'V- V \ 

X. 


IT” 

1 ^3 

Itera- 

f time 


|tion 1 

1 

1 tlon) 

1 tions. 


1 

! I 


lion/ 


I -4 

$ 10 in, 

1 

1 


I 

1 i J 


1 Sq.in. 


i 

1 


10^ 

Initial 

5.000 

3.846 

4.273 

5.485 



1 

800.00 




18 

20 ; 27 



Optimal 

9.090 

11.153 

21.356 

•2,231 





Initial 

5.000 

3.846 

4.808 

7.726 



2 

ao^ 

800.00 





12 

15 { 39 



Optimal 

8.106 

11.38 

14.0 

3.892 





Initial 

4.000 

3.846 

5.495 

16.28 



3 

800.00 





6 

06 ! 50 



Optimal 

7.122 

12,054 

9.000 

6.306 




'4-0% 

Initial 

4.000 

3.846 

6,410 

19.900 



4 

■800.00 





13 

13 ! 2C^ 



Optimal 

6.190 

15.714 

20.218 

9.930 





Initial 

5.000 

3.846 

4.273 

5.485 



5 

10^ 

900.00 





27 

21 ; 03 



OptiDial 

9.238 

23.993 

20,84 

2,208 





Initial 

5.000 

3.846 

4.808 

7.726 



6 

20% 

900,00 





14 

13 : 53 



Optimal 

8,108 

22.28 

17,34 

3.^76 





C'';’'. Initial 

4.000 

3.846 

5.495 

16.280 


, 

7 

30% 

900 .«to 





13 

13 s 23 



Optimal 

7.011 

24.875 

25.615 

6.268 





Initial 

4.000 

3.846 

6.410 

19.900 


15 ! 50 

8 

40% 

900.00 





14 



Optimal 

6.265 

25.0S5 

27.910 

9.820 





Initial 

5.000 

3.846 

4.273 

5,485 



.9 

\0% 

1000,00 



30.62 

2.193 

24 

21 ! 08 



Optimal 

9.204 

30,00 





Initial 

3.000 

3.846 

4.808 

7^726 


12 } 58 

10 

20% 

1000.00 




3.846 

17 


Optimal 

8.109 

30, 200 _ 

28.000 



11 

30% 

Initial 

1000.00 

4.000 

3.846 

5.495 

16.280 

20 

14 { 56 



Optimal 

7.161 

30.35 

31.0 

6,195 




12 40 ^ 1000.00 


Initial 4.000 3.846 ;6.410 19.900 
Optimal 6.276 32.000 41.8 9.747 


18 


13 t 11 



Data 


TJBIE 7 

Initial thickness of the strip being rolled = 0.063” 
Maximum allowable deflection = 5 x in 


inches 


S.No.l $ 

Tzmz 

Design Vector 

1 Cbi.fn. 

X Deflec-i No. 0 : 
1 tion ofj itera- 
‘48 Rolls J tions 

1-4 I 

p lii f 

fjrCorip 

tijBB 
1 !!in 

l 

1 

][Eeduc- 
I tion 
][ 

1 

-1 1 

i i 

l ! 

J L 

^1 

ia. 

i Xg )i x^ ][ (Power 

1 . 1 . 1 Eefjaire 

JTon/sqJlon/sq.J Ton-in/ 
1 tti.^ in f rev. 

o 

Initial 

Optimal 

7.441 

7. *26 2 

7.464 6.95 

1.2 X 4.8 X 

~3 -3 

10 10 

790.32 

374.24 

2,804 

12 

3.005 

1— jftb I'l— 1 mi m 

12 i 

2 20% 

Initial 

Optimal 

7.723 

7.681 

2.082 2.416 

2.96 X 4.44 X 

-3 -3 

10 10 

774.17 

646.88 

3.580 

3 

3.640 

03 : 


Initial 8.045 0.998 1.377 980.64 3.973 

3 30^ 

Optimal 8.013 2.57 x 2*88 x 918.18 4.015 

-3 -3 

10 ^ 10 

8.500 0.794 1.321 1768.81 4.13 

13 14 J 

8.183 4.2 X 1.06 X 1168.73 ' 4.397 

-3 -2 ' 

10 10 


Initial 


Optimal 


7 08 J 


iter 
5 Sec. 

27 

57 

14 

29 
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lf.1 CONCLUSIONS 


Analysis of Cold Rolllri;'; Process 


The present analysis gives reasonable results 
as far as roll forces and torques are concerned. These 
results are expected to come very close to the experi- 
mental or actual values if the modified roll radius, 
after considering the elasticity of rolls, is used in 
the analysis . 

The pressure distribution along the work -roll 
interface agrees veil with the general shape observed 
during experimental investigations 

The deformation pattern and the flow of materia 
gives a fair idea about the elastic and plastic deforma- 
tion taking place in the strip. 

If. 1.2 Design of Cold Eolling Process 

The following conclusions can be drawn from 
the optimum design study of single pass rolling - 

(i) In thinner strips, the tensions are more effec- 
tlire in improving the quality of the product 



while in the case of thicker strips, an increase 
in the roll radius giv® a "better Quality of the 
product . 

(ii) in tlie case of thinner strips, when the total 
power available is the limiting factor, then the 
power supplied through tensions is a better way 
to improve the quality whereas in the case of 
thicker strips, the pov/er supplied should be 
through the torque applied across the rolls for 
improving the quality. 

(iii) In order to reduce the running costs, i.e, to 
reduce the pcjor, the tensions have to be made 
zero. 

In the case of multipass cold rolling, the 
best vray to improve ■ the quality of the product is to 
have S2cme|dimensions for the rolls of all mills, almost 
equal tensions at all the places a.nd to give the maxLnium 
possible reduction in the first pass itself. This 
conclusion is based on the results of one example and 
is to be supported er.tli further numerical results. 



yy 


4.2 EECCOrilEIffiATIOl!' FOR FUTUEE WOEIi 


Tlie finite element analysis can te improved 
lojr tailing; liigher order terms in the displacement model. 
Besides, the deflection of the rolls should he taken 
into account to obtain better results . This analysis 
can bo oxtendod to rectangular bars and rods, uhere 
spread, is considerable, by considering a 3-dimensional 
element. 


Tlie f ini bo element analysis can also be used 
for tlioimal nnn.].ysis to find the temperature distribu- 
tion and J.i'Mico tiic material properties of each element. 
This in J orient ion, in turn, can be used to make an 
elasto-plastic analysis of a hot rolling process. 


In tliG case of design problem, the other types 
of ob.ioctlvo functions like shape function, rate of 
production oto . can cdso te considered. Besides, by 
oonsldorinc the ohanBS in friction due to change in speed 
of rolls , .-elooity can be tahen as a design variable . 

Similarly, the rate of production can be taken 
as the ol.;,joct.LyB function in the case of tandem mill. 
Also the :protslem can be solved for finding the optimum 
design of tandem mills having more than two single pass 

• , T, o in -cries . In the present study , rolling 

rolling mi -Lis in -cries. j:- 

mills having two rolls only are considered. The problem 
can bo extondocl for the design of rolling mills having 
back up ro3.1.c , W-oh are .ulte extensively used in 
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APPEI'iDlX A 


"VERTICAL DiSPLACEi-iEJIT. (Ay) OP A NODE IN CONTACT WITH 
THE ROLLS FOR DISPLACEKENT ( AX) IN X DIRECTION 

Let the coordinates of any node i he y^), 

and the radius of rolls he R, Let the inrut and the 
output thickness of the strip or the "bar rolled he 
and respectively. 

From Fig, (3B) 

= CD + DQ 

= CD 'i- >/r^ - OE^ 


= CD + V R (h^ - h^) (h^ - (A- 1 ' 

The new coordinates after a movement of x are 
■ ^^i +Ay) 

Therefore 

E^P = OP - OE 
/ 

= R - OE 

OE^ r2 ^ b'e ^ 

iff \ 

BE = DQ=CQ- x) 


CQ 


OE 


Hence 


CQ 



From Eq, (A-1), one gets 



(A-2) 



f t 

B E ?^n Eq, 


‘^A-'3) is defined in Eq, (A-2) * 


(A-3) 


This can he implemented on computer to calculate A j 
for a unit incremental displacement 4 x. 



APEEII3DIX B 


DERI7ATI0N OF FONCTION P-. Al© P^ 

3 5 

Inserting the values of s'*‘ and s“ into 
Sq. (3.7) gives 


P = R k 


0 

/ 


h 


h. 


(1 - 


^ ) e 






h- 


(1 - 


n 


- H) 

e ^ d ^ 


As sume 


a = 


Td = 




h. 


1 - 3^/ k 

1 - (ff / k 

1 - h ^ 
h 


1 

For ease of calculation put x 


_ vTH? 

ho 




Hence dx 


h. 


d J25 


One gets from Fig, (1), t|ie thickness of transverse 
slice at distance x, 

ho + R 


h 


h. 


h/- 


= 1 + x" 


3 

k . 



from this one gets 


h 

h 


h 


h 


h^ • ~ 2') (1 x^) 


1 "1 ""2 

-1 


p0l _ 2 a tan** x and = 2 a tan"'' 


/ > 




Now P = R k (1 ;;^) 

k 




1 - 2 a tan 
b (1 + x^) ^ ^- 


-1 




"1 


R “ ^ + J (1 - r) (1 + x^) e 
"n 


2-v 2 a tan"^ x 


\/ ^ 
T 


d X 


or P = 


= k ( 1 - )\/~R Ch- - h„)[ v/T 

k . T , 2 L 


r 


X, 


'n ■ 


-1 


J b (1 4“ x^) e^ ^ ^ dx 


/ 1 - r ^2 a tan 


-1 


+ (1 - r) ' — e‘ 






/ ■ (1 + x^) e" 


^•1 

2\ „-2 a tan x 


dx 


] 




Comparing above equation with Eq. (3*1l) one gets 




(a, r, b) = b’^^^ J b (1 + x^) e 


-1 


2^ 2b. tan X 


dx 


o 


-1 

^ 23. tsil f N H 

+ (1-r)''k2^ e ^ J (1+1®) 

g-2 a tan*’' s _ 



Roll torque function can be derived by inserting 
the value of s*** and s" into 'Eq. ( 3 . 10 ) 






G = R R k ( -I . n )l y h. ^ 


k L - hj I’ ® 
o 


0 


+ ' iT ® 4 J0 


/ 


for ease of calculation a new function z given by 
0 is used. Writing ^ ^ and /U h, in tenos 
of it and. substituting in above equation one gets 


IT 

G = E k ( 1 - 4 ) (h, - hg) I / b y C1 + x2) , 

t o 




■k 


~1 

„2 a tan x , .t* \ 

xe dx+(1~r)e 


2a tan x 


1 


,-1 


/ (1 ''•2'' '' ■'“■2 a tan ' X 


+ x^) X e - “ dx 



Comparing above equation for G with Eq. (3.12) one gets 


f j - (a, r, b) - I hg"* 


{1 + ^2, ^ ^2 a tan-ir^^ 


+ (1 - r) e 
e-2 ^ = ix 


■“i 

2 a tan“ x 


] 


^ / (1 + x^) . 


These values for f^ and f^ can be calculated by numerical 
methods for any values of non dimensional constants 


a, r and b. 



O'OO ooooooooo 


rlM=C 

;‘,.L;'a=L)f;LX 

1C C’JI.TINUE 
Pk!NT2,,Nl 
C/LL I'JCDIS 

form stiffness matrix 

CALL FORMK 

scLvt The equations 

Call soLvt 

CALCULATE THE STRESSES IN EACH ELcMcNT 

CALL STRCSS 

CALL SIG8AR (ANETAiMMfNIf Nil } 

IFiMMoLQe i}GCT02C 
DC LX = DC LX 1* I ANC T A- 1 .. 0 ) 

GUT01C2 
2C CGNTINUd 

! F (D"L'%GT«C .CDGOTLIOI 
'->':nnic2 

101 "■'.L';=CcOi 
IC2 CJMIKUc 
■!I=NI + 1 

IFlNUGToAPR )GLT0;>G 
CulQlC 

: c CiiMTiNur 

1 FiiRMATC 91.;) 

2 PORMAT(22h NUMBER OF I TERATI CNS=i lE, 20X, 3HMM=, 12 ) 

CTCP 

• ND 

5UBRCUTIN „ GCATA 

CiJMMCil/CONTR /TITL: ( 12) f NP,NE t NDF, NCNt NLD, NMATt MS ZF ,NPREf LMCN 
CriMMON/CONTS/CORDtSOOf 2) ,N0P(472t3), ORT( 1, 2) 

COPMCN/CONTT/IMP (600) ,Xt600) |B(600) 

Ci:(NMCN/C0NTU/DELX,RAC,Hl,H2,AD,Y 

C<JMMCN/C0NTV/NH,NV 

*■ 

READ AND PRINT TITLE AND CONTROL 

Rt AD7,TlTLf 
PRINTIOO, TITLE 

R •' AD 1 f HP I Nf » N LD f NOF I NMAT f N PRE 
P R IN T 1 1 NP f Nl; f NLD fNDF j NHAT j NPRE 

READ AND PRINT MATERIAL DA3A 



c 

c 

c 


c 

c 

c 


R ^D.3,(N,{0RT{N,I},I = 1,2),l=1,NMATJ 

^’Kinioe 

PRINTS, ff^f(CRT{N, n, 1 = 1,2), N=1,NMAT} 

GCNSRaT^ nodal pci NT DATA 


k.,adi,mv,nh 
R‘ AD€,STi,:PH,STt*PV 
■*.2=N V-l 
DD£OI=l ,NV 

1 — O c c 

Du50J=I ,NP,N V 
C URD { J, 1 ) =A 1 ^STt PH 
CORDC J, 2) =A2*STfcPV 
i', 1 =A 1 + 1*0 
’- C CONTI NUl 

C2=A2“io 0 

6C CUNTiNUe 

G-NcPAT.. fcLLMENT DATA 
P.ADl,.Tn,II2,n.,JJl,JJ2,JJ;^,il,I2 

■iHH = NH-l 
'V\/=NV“i 
DUf?OI=l ,-!HH 
: Jl=X!l 
!: J 2=!!2 
! 

Jl 1=JJ1 
Ji 2=J J2 
\i ^ vJ D 
D:J?0J=1,NVV 
'lup{ II, 1)=IJ1 
'0P{U,2)=I J2 
MUPd 

;!UP{I2,1)=JI1 
'1()P( 12,2) =JI2 
>ap( 12,,;.) =J13 

r 1=1 1+1 

I 2=1 2 + 1 
7 J1 = IJ1 + 1 
! J2=IJ2+I 
’ J*! = IJC+1 
Jl l = JIi + l 
JI2=JI2+1 
JI2=JIX+1 
TO CONTINUE 
I 1=1 1+NVV 
! 2=1 2+NVV 
I11=II1+NV 



o o n o o o 


J 2 =II 2 + IMV 

■J Jl= J Jl + NV 
J J 2 =JJ 2 + N'V 
I J 2 = J Ji + Nv 
5C CONTI NUr 

'•H V=NHH#NVV *2 


RtAD AND PRINT BOUNDRY DATA 


PRINT 104 
NSZF=NP*NDF 

RJA04,(INP(I},I=1,nszf> 

R : A D 1 0 , Dfr L , R AD f H 1 , ri 2 , AD 
Re AD 11, Y 

PRINT10,D£LXfRAD,Hl,H2fAD 

PRINTlltY 


PRINT INPUT DATA 


Pf<lNT102 

PRINT2, {N,(CGRD{t,,M),A=l,21,N=l,NP) 

PRINT (N, (KCP(N,P) 5P'=l,3)tN=l,N£) 

;CC CONTI^U;., 

1 FORMA T( 91 ; ) 

2 F(iRMAT(3{IS,2F10e4) J 

0 F u R M T ( 2 '> 1 ) 

format (,. 01 1) 

5 FORMAT! FGo A) 

T: pOPMaT! 2F12® 8) 

7 FORMAT! 12A6) 

8 FORHAT(iT,2£l0.3) 

9 FORMAT! lyjGlCcS) 

1C FjRMAT! nF6«4) 

11 FORMAT! . -I Co2 ) 

ICO FORMAT! IH ,12A6) 

1C2 FORMAT! I! H NCDAL POINTS) 

1C3 FORMAT! 9H oL MENlS) 

104 FORMAT! 20H BCUNDARY CONDITIONS) 

ION FORMAT! /3CH GIVlN COMPONENTS OF VECTOR X) 

1C6 FORMAT! /?CH GIVFN COMPONENTS OF VECTOR 8) 

108 FrjRMAT{2CH MATERIAL PROPERTIES) 

RPTURN 

■ oo 

SUBROUTINt: INCDIS 

COMMGN/CONTS/C0RD{300,2)tNGP!472,3)fORTI 1,2) 

CDMMON/CONTR/TITLi' (12),NP,NE,NDF,NCN,NLD,NMAT,NSZF,NPRe,LMCN 

COMMON/CONTT/IMP !600 ) ,X{&00) ,8(600 ) 

C0MH0N/C0NTU/DELZ:,RAD,H1,H2,AD,Y 



: '^'MC^.!/CaNTY/NH,NV 
.■■i' = l,NSZF 

! F (li^FCl )c'fC,>0)GOT010 
: ( I) =Co 0 
3 il ] 0 1 ':i 
1C 3 ( i j =Co 0 
r. CONTI NUi; 
n;t20I=l,NP 

2C C(JRD{I,1)=CGR0{I ,l)+D‘gLX 
do‘;'OI=i,np,xv 

I F (CORD (I fl)„L£e AD)GCTC30 
ai=SQRT(RAD=^(Hl-H2)-(Hl-H2)**2/fc,0) 

:i=AC+ABSCXl 1) 

I F (CORD? I , 1) sjGTo XI )G0T025 
Y1=RAD+H2/2«C 

\P=SGRT(RAD*^2-(C0RD{I, 1)-X1J=»#2) 

YN=Y1-ABS(AP] 

Y{2*I}=YN-CCRD(I,2) 

IMP( 2*1 J = 1 
;; ( 2*1-1 )=C«C 
JMP(2*I-1)=1 
RRi'lT6C,I 
Girro:’''c 
2 ;- a.tM'lNU. 

P. {2*I) = 0., C 
I HP { 2*1 )=C 
R(2*!-l)=CoC 
1 MP( 2*i “1 ) =C 
pRirn^Oti 
;;C corofiNuc 

^-C FORHATCI ">,*TH fiQD*;' HAS COHSQUT UF THE ROLLS*! 

PC FORMATC* iyp&s=*,iooi 1 ) 

6C FURMAK iCf^TF NODI- IS BfcTwEEh THi* ROLLS*) 

ROTURN 

iiC 

RUBRCLTiNfc FDRMK 

CQPMON/CONTR/TITLF J 12) ,NPtNE»NDF»NCM,NLD,NMAT,NSZF,NPREtLMCN 
C JPMCN/COKTS/COROOOO,2) ,NCP(<f72,3),ORT( 1,2) 

C JPMCIn/a8G/SK(510, 14) 

C JPMCN/ABC/L-STIFK{6,6) 

CaPMON/BCD/PB,NNN 
DI PENSION KK{G) 

DDC0CN=1,NSZF 
DD30CY=1,PB 
jOC 5K{M,M)=CoC 
HBI = C 

DU5C0N=1,Ne 
CALL STIFT2(N) 

KK{2)=2*NCP(N,1) 

KK{4)=2*NGP{N,2) 
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KK(i) = 2 ^riCF(,N,?| 

KKU ) =KK(2)-1 
KK(.. }=KK( 4)“1 
KK{ ;) =KK( 6)-1 

ni.>} 0 ci = i ,6 

1 i =KK(I ) 
n'l 4 CCJ=l ,6 

I F {KK( J), LT* JI )GCTOvOO 
J J=KK(J)-II +1 

SK(II,JJ)=SK(II, JJ)+cSTIFM(If J) 

! F CMBIr- GT« J J JGOTG ^-00 
MBI= JJ 

-iOC CONTINUE 
50 C CQNTINU:'.* 

M 3 =MBI 

FRlNT 600 fMBI 

oCC FORMAK* BANDWIDTH==«^,I 5 ) 

RETURN 

FND 

SUBROUTiMt STIFT 2 (NJ 
CC 1 RM 0 N/A 8 C/E 5 TIFF( 6 , 6 ) 

C OFHC^'/ aBC/ A (I. , 6 ) , B ( 3 , o ) 

:., LL ST 8 AR(NJ 

ESTIFH IS STIFFNESS MATRIX 

D 021 CI= 1 ,C 
0021 CJ= 1,0 
■ Cl'iFN{l,J)=C «0 
')a 2 lCK=l 

21 C oTIFMC I , J) =,.:STiFM{i,J) + 8 {K,I)*A(K,J)*AR£A 

R .TURN 

ilC 

‘^UBRnuTINi: STBAM{N) 

CU.WMQN/CQNTB/TITLF C 12 } , NP, NE.NDFtNC^ ,NLD, NMAT, NSZF,NPRE,LMCN 
C 0 RMCN/C 0 NTS/C 0 RD{ 3 C 0 , 2 ),NCP{ 472 , 3 ),QRT( 1 , 2 } 
CJMM0N/C0NTU/DELK,RaD,H1,H2,AD,Y 
Crjf^MCN/ABD/A (3 ,6 } , B ( 3 , 6 ) 

COMMCN/ABF/ARFA 

C'JRMnN/ABF/FCRCi ( 410 , 4 } ,STRAN( 410 , 3 ) , S IGB( 410 } , INP ( 4 10 », cFSTC 410 } 

Di MANSION USTIFM{ 6 ,»j) 


DCTCRMlNc 6LEMENT CONNeCTION 


I =H0P(N,1) 
J=N0P(N,2} 
K=N0P{H,3} 
L = 1 



o o o o o o 


c 

c 


SET UP LOCAL COORDlNArC SYSTEH 


' J=C CRD (J,li -CORD (1,1) 

^.K=C CROCK,!) -CORD (1, 13 
BJ=CGRD( J,2)»-CORD(l,23 
?K=CGRD(K,23-CORO(I ,23 
K B A = ( A J* BK- AK=i'B J 3 / 2, 0 
■■ R.vA=/;BS{AR^-A) 

I F (ARr-AeLL-« CoO)GCT0220 

FORM STRAIN DISPLACLMFNT MATRIX 

k I 1, 1) = BJ-BK 
P. { 1, 2 3=0* C 
A(1,23=BK 
A ( 1,4)=0»C 
A{ 1, *'3=-BJ 
Atl,€)=OeC 
A{2,1)=0„C 
k ( 2, 2)=AK-AJ 
A{2,33=CaC 
•'U2,4)=-AK 
(2, ■)=0,. C 
A (2, J 
, 1)=.,K-AJ 
, 2) = BJ“BK 
■ ( . , ' ) =-AK 
'■ (3,B)=BK 
A (2 , C)=AJ 

002011 = 1 ,.. 

01,1201 J = 1,6 

2C1 * (I, J)=a(I,J)/(2,0=<'ARCA) 

FORM STRSSS-SRRAIN MATRIX 

' FClNP(N)ebC«l)GCT0202 

CaPM=ORT(L,l )/( { U 0+0RT(L,2) 3 * ( O-ORTC L, 2 3 0 3 3 
STI FM{ l,13=C0MM=i‘{U 0-CRT(L,23 3 
CTIFMt l,2)=CaMM=»ORT(L,23 
STIFMC 1,3 3=C«0 
' STIFM(2,13= = ;STIFF( 1,23 
"STIFM{2,23=:STIFM{i,l3 
CTI FP{2,33 =U0 
' STIFM(3, 13=0.0 
: S7IFM(3,2)=C. 0 

STiFM{3,3 3=GRT( L, l)/( 2.0*(l.0+0RT tL,2) 3 3 

GOTO 204 

202 :uM=ORT{L,l)/(leO+ORTtL,233 

C0M1 = (1.0-0RT(L,23 3 / ( 1. 0-2.04QRr ( L,2 3 3 



c 

c 

c 


CUF2 = 0RT( Lt2)/(UO-2cO*ORTCLf2) ) 

ilP=OeO 

n l*0'^^^T{L,2))*HP/(ORT(L,i)*3.0n*SIGB{NJ*#2/3. 

0 J = FGRL,.{ )+F0RC.r;(N»2J + F0RCi.'(N»4) 

') i r G': =F ORCE ( N , U -D J/?o 0 
n.'iIGY=FaRCHN,2)“0j/3o0 
STI FM ( I , U =€0^* (CCFl-DSIGX*=f2/CQMS) 

STI FM( i»2) =C0M^{C0R2-0SIGX«0S1GY/C3I«S) 

STIFM{ 1,2) =C0M*{-D5IGX*FCRCS=(Nt3 J/COMS) 

STIFM(2,l)=.STIFM(i,2) 

STI FM( 2,2) =COM* (C0Rl-DSIGY*#2/C0MS) 

' STI FM( 2,3 J =C0M#(“DSIGY*F0RCfc;(N,3 )/COMS) 

: STIFW{; ,1)=.;STIFM{ 1,3) 

• STIFRC3,2)=cSTIFR(2,3) 

' STIFR(3,3)=C0M#(0*, -FaRC£{N, 3) **2/0085) 

tt1 MINUS 

D0205I=1,3 

D0205 J=l,6 

B (I, J)=0, 0 

D0205K=1,3 

BCI, J) = B{2, J)+ESTIFRII,K)*ACK, J) 

G-;"TURN 


ERROR tXiT FOR BAD CQNNECnON 

22C PRI’ITlCCjN 

ICC ^:jRNAT{;.3H Z RC QR : EGaTIVE ARtA ELEMENT N0,14/21H ‘XcCUTION TtRMi 
It '1 D) 

STOP 

-HD 

SUBROUTIMA SCLV- 
DiMtNSiCN XX (6001,Ii ntOO) 

C,lMMON/CONTB/TITl^' (12) ,NP,NE,NDF,NCN,NLD,NMAT,NSZF,MPRE,LMCN 
CLlYMaN/C0NTT/lMP{600) ,X(600) ,B(600) 

CUVMCR'/rtBG/SKCGiO,!^) 

CUNMON/BCD/MBtNNN 
CUMMCN/XYZ/BB (600 
n;jcON = l,NSZF 
IF(IMP(N)6cG<il)GCTOaO 
nu-<iOL = l,MB 

I F(SK(N,L)oSG.0,0)GCTC40 
I =N+L“1 

' F{IMP{n.feG«OyGGTOAO 
B(N)=8(N}-SK(N,L)*X(n 
4C COMINUe 
nO'50L=2,MB 
I =N+1-L 

IF{I®LSo0)GCT050 
TF(IMP(I)«HGi.O)GCTDiO 
B(N)=B(N)-SK(I,L)*X(I ) 



C CnhTlNU?^ 

*:.C Cj-NTINU:-; 

C FGRkA 5^D RiDUCTlON CF MATRIX (GAUSS ELIMINATION) 

!'iil79CN=l,NSZF 
" F (I MP( N)« -GoDGCTOyqO 
'■»):i78CL=2,MB 

' F CSK(N,L)*FC*0«0)GUT0780 

T =N+L-1 

.TF(IMP{I)»teCal)GCTaV80 
r. = SK(N,L)/SK(N,lJ 
J=C 

DQ75CK=L»MB 
J = J+1 
I I = I + J”* 1 

J F(IMP{ II ).E(;.1)G0T0750 
SK(I , J) =SKCI t J)-C=S=SK(N,KJ 
15C CONTINUE 
SK(Nf L) =C 
780 CONTINUE 
79C CONTINUE 

C FORWARD REDUCTION OF CONSTANTS (GAUSS tLlMlNAIION) 

D08'iCN=l,NSZF 
; F(lMP(N)e,;.G->l)GCTQ:.70 
n.J32CL = 2,M8 

;> F { S K { N , L ) • . Co Oo C ) G CT08 20 
; =f\+L-l 

TF(lMP{I)o:..C:al)GCT0ci20 
5(1 )=BtI)“3K (N,L)*B(N) 

L2C continue 

B(N}=B{N)/SK{N,15 

e: C CONTINUE 

C SOLVE FOR UNKNOWNS BY BACK SUBSTITUriON . 

nO£6CM=2,NSZF 
■■' = NS2F+1-M 

■:F(IMP(N)oECo1)GCT0860 
00&‘:.CL = 2»MB 

: F {3K{N,L)o;iGoG«C}GCTQ850 
K. = N+L-1 

J F (I MP( K)oKC;ol)GCTQ650 
3 (N) =3{N)~SK (Nf L}*BtK) 

C.:'C CONTINUE- 
hhC CONTI NU.' 

1? F0RMAT(/07H COMPUTED COMPONENTS OF VECTOR X) 

JJ=C 

90141=1, NSZF 
F (I MP( n*CG« 1 )GCT014 
):ii)=8(i) 

J j=j j+i 
EX{ jj)=x(n 



o o 


FST(I) = ‘FST(D+YZX 

,..cc : 

F .Jf<C (1 ,'?)=FCRCl. {I,‘.) + FINT(I,4} 

P-IRC (i, J)=FCRC*,(I, J) + FINT{I,J) 

.''.TRANli. , J)=STRAN(I, J)+SINT{I ,J) 

1-,.C CufvTlNUt,, 

PRINT DiSPLACcMrNT OF NODES 

Dcr :ci=i,KP 

C U Ri) ( 1 , 1 ) =CC RO { ! , 1 ) +DI S { 1 , 1 ) 

CiiF;D(I,2) =CCRD{1 t2)+DISI2,I) 

120 CONTINUE 
PCTURN 
I i C 

CUBRCUTiNi SiGBAR(Ah5tTA,MH,Nl,NlI} 

CDFMCN'/CaNTR/TITLI (12) ,NP,N£tNDF,NCN,NLD,NMAT,NSZFtNPRE,LMCN 
CnFMCN/CONTS/CORDOCOf 2)f NCP(%72,3)»ORT{ 1*2) 
COFMCN/CaNTU/Db“LX,RMD,Hl,H2,AO,Y 

Ci >FM0^4/A8F/F0RCr (AlO »A-) »STRAN(<i-10,3 J tSIGB( 410 ) i INP (410 J t ; FST( 410» 
COFMON/CDE/DIS(2,'00) 

CUFMCN/BCD/FBtNNN 
D01CN=lf NE 

A=( FORCE {N,l )-F0RCE(N»2) }=i'=«'2 
P = {F0RC ,(N»2)-FURCE{N»4) )**2 
C = {FCRC:.:{Nt4 )-FORCE(N, 1) ) **2 
1C SIGB(N) =SGRT { (A+B+C)/2*0+3«0*F0RCfe(Nf 3}**2) 

11 ClJNTlNUO 
SiG8F=0o0 
I =C 

D.J2CK = l,Ne 

r F{INP{N)*uG«l)GGT020 
J F (S1G6(N). LO 0 SIGBFJGQTO 2 O 
Si.GBF = SIGB{N) 

I =N 

20 CONTINUE 

I F (MFeL:G.C)GCTQ4C 
I F (SIGBF.« G'£« C.995*YJG0T03O 
GOTO.‘;0 
OC !liP(I) = l 
LMCN=1 
PRINTTOf I 
GOTO 11 
40 CONTINUE 

I F ( SIGBMoGE- C»9g5*YJGCT050 
AN5TA=Y/SIG8F 
GOT060 
5C INP( I )=1 



LKCN=L 
>1 M = 1 

Gmnu 

,;'C .INTIN U : 

7C =URMAT C 15,2-;H TH HAS YluLDgD) 

1 p (Wlo QoNPR;:iG0Ta5C0 
F (WI® ■ "Q*MI )G0TG'; 0G 

G iJ T 0 V C O 
70C CUN'flMjG 

'U i = Nr I +NKN 
PRIWTlOl 
G 

C WRITE ALL STRESS CCHPONENTS 

C 

ni.!f,ocN=i,hs 

PKliHl 1 1,A, (FORC*. (N,I ),I = l,4)fSIGB(N), (STRANIN,! ),I=l»3»i£FSTtM» 

OOC CCMlNUfc 

ICl FORMAT ( lXt*cLGMi;NT=J=,'iXi*X-STPPSi=i',&X,*Y-*STRESS*,5Xt*XY-STRESS*,6X, 
L*Z“STRi.:SS ‘3FFrCT4STR::SS=«‘,6X,*X-STRAIN*,6X,*Y-STRAlN*f5Xf<=XY-STRAlM 

2 FFF'CTo STRAIN*) 

111 FORMAT! 15, iiFl4«4,4F14«6) 

PRINTICO 
00120^ = 1, NP, 2 

PRINTl 10,M, (0IS( J,MJ ,CCRD{ M, J) , J= 1,2 ) , MLt C DISC J, ML ) f CORD! MLt J ) f J = 1 
1,2) 

12C CONTINUE’ 

IOC FORMAT! ///,lGX,29HDiSPLACeMENTS AND 'COORDINATES ) 
lie FORMAT! 2{ 15, 4Flia6) ) 

LMCM=0 

70C CONTI NU; V — 

R A TURN 

■ ND 


\ 



o o o 


'^is(2,j)=){(in+i) 

9C j;i = ni+2 


Calculate titHENT forc6S 


n'ji3ci=i,Nf 

FINT{I,4)=FCRCSf{i,4) 

FxNTJI, j)=FCRCc{i,J) 

SiNTd, J3=STRAN(I, J) 

UC CONTiMUSi 

1'>02CONC = 1 ,NS 
CALL STBAP(P«C} 

D026CI=1,NCN 

M=NOP(MC,I} 

X F{M*liQ®0)GCT0260 
K = {I-U*NDF 
DU24CJ=1,NDF 
IJ = J + K 

24C R {1J)=0IS{ J,P) 

26C CONTI NUfc 
i A=K+NDF 

nn,:;oc! = i, 2 . 

FORCv (NC,I )=C„0 
ftr;?' (oc, I ) =CcO 
J Ll .. C C J = 1 j I 

F JFC ( !C, I) = FaRC .{NC,I )+BB(I ,J}=4=R{ JJ 
STRAN( WC,U =5TRAN(NC,n+AA(r J) 

:;CC COMTIfJUc 

2CC continue 

‘ FdNPa )<,LLoO)GCTDi 20 
CnP=CRT( 1 ,1) /{1«C+0RT{ 1,2) ) 
rrjM = CRT( 1,2)/ ( 1,0-2* 0*0RT( 1,2) ) 

HF=0,C 

Cl)P5 = 2o 0={=(lo C+2.0*{ UO + GRT 11,2)) *HP/(ORT( 1,1) *3.0)) *SIGB{i^i) **2/3,0 
nj=FINT{I ,1) +FINT(I ,2)+FlNT{ 1,4) 

DSIGX = FlNT(I,l)-DJ/..,0 
D.':1GY=FINT{I ,2)-DJ/: ,0 
0SIGZ=FINT( I ,4)-DJ/A,0 
CAB=CnM*(CCf'l-DSIGK*DSIGZ/COFS) 

CAC=CaM*{COP 1-DS1GY*DSI GZ/COFS ) 

CaU=C 0M*(-0S1GZ*F1NT{1 ,3)/C0f^S) 

FORCEli ,4 )=CaB*STRAN{I,1)+CAC*STRAN{ 1,2)+CAD*STRAN(I,3) 

GOTOnSO 

i2C FORCrd ,4)=CRT( 1,2)*{FQRCE{ I ,l)+FORCt-(I, 2) ) 

230 CONTINUE 

005001=1, NE 

Z YX=STRAN( I ,l)**2+STRANtI ,2 ) **2+0, 5*STRANt I, 3) **2 
YZX=SQRT(2,C*ZYX/3.0) 





ono ono ooooooooooooooooo 


APP'ND1X«D 


■!P=!iCa UF NCDaL POIMS 
=NC» OF kLt.MENTS 

''B=N 0 o OF RI:STRAINl-D BCUNDRY NODcS 
■! 0 F=Ka 6 OF DEGREES OF FREEDOM Pc.R NODE 

MCr,’=NA/U OF nod; s per element 

^*LL:=NaE OF LCAO CASiS 
■' 1 MAT = N 0 « CF material TYPES 
!SZF=NO« CF uQATinNS IN THE SYSTEM 
L 1 =lCAD CAS" COUNTER 
COORC=NOOAL POINT CCCRDINATE ARRAY 
10 P=~LEMENT CONNcCTiGN ARRAY 
IriAT = r:LL”M£NT MATERIAL TYPE ARRAY 

:)rt=::l*„;hent type matlrial prcplrties 

NBC=Rr:STRAlMD BCUNDRY MODE NUMBERS 
•?FIX = BOUHDRY CONDITION TUPE 
R 1 =LCAD VECTOR 

SK=Rf:CATANGLLAR MATRIX FOR SCUATlONS 
COMMON/CONTR/TITLF i 12 ) , NP, NE, NDF, NCN ,NLD, NMAT, NS ZFf NPRE t LMCN 
C 0 MM 0 N/C 0 NTS/C 0 RD{ 300 , 2 )fNCP(% 72 f 3 J»ORTl If 2 ) 

CfJMMON/CONTT/iMP (feOC) ,X( 600 ) ,B{^ 00 ) 
CUMMCN/CONTU/DEL>:,R..D,Hl,H 2 ,AO,Y 

COMMCM/aBF/FCRCc {^ 10 ,%) , ST RAN( 410 , 5 ) ,S IGB( 410 ) , INPI A 10 J , FST( 4 10 ) 
C 0 MMrN/;;BG/SK( 31 C,l%) 
r UMMON/BCD/MBtNNL 
:. 3 MMCrj/y;YZ/BB( 60 C) 

M 8 = l^. 

■!CN=:i 
LMCN=C 
M n = 1 

Ri,ADl,NNN 

R-aD1,NPRCB 

LCGP ON NO. OF PROBLEMS 

Dr) 4 CCNPR=l ,NFR 0 B 

RfcAD INPUT GEOMfeTRY AND PROP. AND B.C.iS 

CALL GDATA 
-sCC CONTINUE 

CALL FLUNCC. 2 C 00 ) 

DUn=l,NS 
-FST{I)= 0 . C 
IMP (I )=0 
D 05 J= 1»3 

FORCCd, J)=CoC 
STRANd » J)=C 0 O 
5 CONTI NUC 
Nl = l 



